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Seventeen  simply  supported  reinforced  concrete  beams  were  tested  to  collapse 
under  combined  flexural,  axial,  and  shear  forces.  Each  beam  had  a 12.5-foot 
span,  a 9-by-15-inch  cross  section,  and  an  effective  depth  of  12.5  inches. 
Reinforcement  consisted  of  No.  2 stirrups,  6 inches  on  center.  The  beams  were 
loaded  laterally  through  a symmetrical  two-point  loading  system  and  axially 
through  the  plastic  centroid.  Load  was  applied  by  a single  hydraulic  system 
designed  to  provide  a constant  ratio  between  axial  and  lateral  loads  for  the 
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20.  ABSTRACT  (Continued) 

'duration  of  the  tests.  The  two  test  parameters  were  axial-to-lateral-load 
ratio  and  shear-span-to-beam-depth  ratio.  Electrical  measurements  of  beam 
behavior  included  steel  strain  on  the  longitudinal  rebar,  concrete  strain, 
vertical  deflections  along  the  length  of  the  beam,  end  rotations,  and  lateral 
and  axial  loads.  In  addition,  a photoelastic  coating  sheet  was  bonded  to  one 
side  of  the  beams  and  overlaid  with  a sheet  of  Polaroid  film.  The  experimental 
results  from  the  beam  tests  were  compared  with  data  calculated  with  an 
analytical  behavioral  model  developed  as  part  of  this  effort.  The  general 
beam  behavior  calculated  from  the  analytical  model  agreed  well  with  the 
measured  results,  especially  in  the  region  up  to  maximum  load. 
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SECTION  I 
INTRODUCTION 
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BACKGROUND 

Although  a tremendous  amount  of  research  has  been  directed  toward  an  under- 
standing of  the  behavior  of  reinforced  concrete  structures,  there  are  some 
areas  in  which  knowledge  is  far  from  complete  or  is  inadequate  for  some  ap- 
plications. For  example,  the  behavior  of  a reinforced  concrete  beam  loaded 
to  collapse  under  the  influence  of  combined  flexural,  shear,  and  axial  loads 
is  not  completely  understood. 

Early  endeavors  in  the  analysis  of  reinforced  concrete  behavior  were  based  on 
the  theory  of  elasticity  and  the  applicability  of  Hooke's  Law.  This  approach 
allows  elementary  strength  of  materials  procedures  and  an  extension  of  elastic 
theory  to  be  used  for  the  analysis  and  design  of  reinforced  concrete  struc- 
tures. Early  investigators  recognized  that  concrete  is  not  an  elastic  mate- 
rial; however,  the  simplicity  of  the  straight-line  method  was  appealing  from 
a computational  viewpoint.  As  better  understanding  of  the  actual  behavior  of 
reinforced  concrete  structures  developed,  revisions  and  additions  were  made 
to  accepted  design  procedures.  This  evolution  has  led  to  the  widespread  use 
of  the  ultimate  strength  design  method  or  atvenjth  "icthod  used  in  the  1971 
American  Concrete  Institute  (ACI)  Building  Code. 

During  the  past  25  years,  extensive  research  has  been  conducted  on  individual 
aspects  of  reinforced  concrete  behavior  (e.g.,  flexure,  axial,  shear,  torsion, 
and  bond).  Also,  many  investigations  concerning  the  interaction  between  vari- 
ous forces  acting  on  reinforced  concrete  members  (moment-shear,  shear-torsion , 
and  moment-axial)  have  been  conducted.  However,  for  a general  loading  condi- 
tion, concrete  member  behavior  depends  on  the  interaction  of  all  these  behav- 
ioral aspects,  and  there  exists  no  general  theory  or  general  behavioral  model 
that  accounts  for  the  combined  effects  of  all  these  behavioral  aspects. 
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The  state-of-the-art  of  reinforced  concrete  design  is  at  a point  where  a great 
deal  of  confidence  exists  in  the  presently  accepted  strength  design  procedures 
when  well-defined  gravity  loads  are  considered.  However,  there  exist  situa- 
tions in  which  knowing  just  the  strength  of  a member  is  not  enough.  To  accur- 
ately predict  the  complete  behavior  of  a multistory  building  to  blast  or  earth 
quake  loads  or  even  the  ultimate  strength  of  a complex  shell  structure  to  stat 
ic  loads,  the  complete  response  of  reinforced  concrete  members  or  elements  to 
collapse  must  be  known. 

In  recent  years  the  complete  behavior  of  reinforced  concrete  structures  and 
structural  systems  has  been  investigated  from  two  directions:  scale  models 

and  full-sized  prototype  structures  have  been  tested  to  provide  experimental 
evidence  of  behavior,  and  analytical  investigations  have  been  conducted  to 
establish  methods  for  reasonable  prediction  of  the  response  of  complex  struc- 
tures to  various  types  of  loads. 

During  the  past  decade,  there  has  been  considerable  success  in  analyzing  the 
response  of  complex  structures  by  the  finite-element  method,  in  which  a con- 
tinuum structure  is  modeled  by  an  assemblage  of  discrete  elements  connected  by 
selected  points  called  ywder. . The  accuracy  of  the  results  depends  on,  among 
other  things,  the  ability  to  accurately  define  the  behavior  of  each  element. 

If  any  analytical  method  is  to  predict  the  complete  behavior  or  response  of 
a complex  structure,  the  complete  behavior  of  reinforced  concrete  under  a gen- 
eral loading  condition  must  be  understood. 

OBJECTIVE 

In  an  attempt  to  add  to  the  basic  understanding  of  the  behavior  of  reinforced 
concrete  members  loaded  to  collapse,  this  investigation  was  undertaken.  The 
specific  objective  of  this  research  was  to  investigate  the  static  behavior  of 
reinforced  concrete  beams  subjected  to  combined  flexural,  axial,  and  shear 
forces  by  developing  a general  model  for  beam  behavior  and  comparing  the  re- 
sults to  those  obtained  from  an  experimental  program.  Of  particular  interest 
was  member  behavior  during  the  large  deflections  which  occur  after  maximum 
load. 
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The  investigation  included  both  an  experimental  and  an  analytical  phase.  The 
experimental  phase  of  the  investigation  consisted  of  static  testing,  to  col- 
lapse, 17  hinge-supported  reinforced  concrete  beams.  All  beams  had  the  same 
span  length,  cross-sectional  geometry  and  properties,  and  nominal  material 
properties.  The  beams  were  loaded  laterally  through  a symmetrical  two-point 
loading  system  and  axially  through  the  plastic  centroid,  as  defined  by  the 
1963  ACI  Building  Code.  Load  was  applied  through  a single  hydraulic  system 
designed  to  provide  a constant  ratio  between  axial  and  lateral  loads  for  the 
duration  of  the  test.  The  two  test  parameters  were  axial-to-lateral-load 
ratio  and  shear-span-to-beam-depth  ratio.  Three  shear-span-to-beam-depth 
ratios  (3,  4,  and  5)  were  considered.  For  each  shear-span-to-beam-depth 
ratio,  three  axial-to-lateral-load  ratios  (3,  2,  and  0)  were  used. 

The  experimental  results  provided  data  to  describe  the  strength  characteris- 
tics of  the  beams,  including  crack  development,  failure  modes,  strength  in- 
terrelationships, and  hinge  performance.  Also  load-deformation  behavior,  in- 
cluding load-deflection  characteri sties , yield  and  collapse  deflections,  and 
load-strain  behavior  at  various  sections  was  determined. 

In  the  analytical  phase,  a reinforced  concrete  behavioral  model,  which  inclu- 
ded beam  response  beyond  the  maximum  load-carrying  capacity,  was  formulated. 
The  adequacy  of  the  model  was  determined  by  comparing  the  model  predictions 
with  the  experimental  test  data. 
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SECTION  II 
HISTORICAL  Rb'vIEW 


Most  of  the  uncertainties  in  predicting  the  flexural  behavior  of  reinforced 
concrete  members  in  the  range  of  their  maximum  capacity  result  from  the  in- 
elastic nature  of  the  two  materials  involved.  Accurate  prediction  of  the 
flexural  response  depends  on  the  ability  to  mathematically  describe  the  in- 
elastic characteristics  of  both  the  reinforcing  steel  and  the  concrete. 

The  reinforcing  steel  most  commonly  used  has  a distinct  yield  point  and  a 
flat  plastic  region;  therefore,  its  behavior  can  be  adequately  expressed  by 
a trapezoidal  stress-strai n curve.  The  stress-strain  relationship  for  con- 
crete, however,  is  more  difficult  to  determine  and  also  more  difficult  to  ex- 
press mathematically.  Discussions  of  the  stress-strain  relationship  of  con- 
crete and  the  stress  distribution  in  the  compression  zone  of  flexural  members 
have  appeared  in  the  literature  since  the  late  1800s.  Early  design  proce- 
dures for  reinforced  concrete  were  based  on  agreement  between  calculated  and 
experimentally  determined  capacities  (the  same  method  the  presently  accepted 
ultimate  strength  design  method  is  based  on).  Therefore,  studies  of  stress 
distribution  in  the  compression  zone  of  reinforced  concrete  members  have  par- 
alleled the  development  of  ultimate  strength  design  theories.  Investigation 
of  inelastic  concrete  stress  distribution  has  been  approached  from  two  direc- 
tions. One  approach  has  been  to  determine  the  concrete  stress  distribution 
by  analysis  of  the  observed  behavior  and  the  ultimate  strength  from  tests  of 
reinforced  concrete  beams  and  columns;  the  other  has  been  to  determine  the 
stress  distribution  by  direct  measurement  of  strain  on  plain  concrete  speci- 
mens. 

One  of  the  first  rational  methods  for  designing  reinforced  concrete  members 
was  reported  by  Koenen  (ref.  1)  in  1886.  It  concerned  the  analysis  of  simple 
reinforced  concrete  slabs  subjected  to  bending.  He  assumed  that  there  was  a 
straight-line  distribution  of  concrete  stress  in  the  compression  zone,  the 
neutral  axis  was  at  middepth  of  the  section,  and  the  concrete  carried  no  ten- 
sion. Koenen' s work  was  followed  by  the  development  of  many  more  flexural 
analysis  theories  on  reinforced  concrete.  Hognestad  (ref.  2)  presented  the 
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highlights  of  several  of  these  theories  published  prior  to  1951.  Figure 
1 presents  the  basic  assumptions  used  in  the  theories  discussed  in  this 
section. 

Among  the  early  theories  published  on  flexure  (following  Koenen's)  were 
other  straight-line  method's  by  Neumann  (ref.  3)  in  1890,  Coignet  and 
Tedesco  (ref.  4)  in  1894,  Johnson  (ref.  5)  in  1895,  and  Ostenfeld  (ref. 

6)  in  1902.  (See  figure  1.)  In  the  Neumann  and  Ostenfeld  Methods,  the 
concrete  was  permitted  to  carry  some  tension.  The  Coignet/Tedesco 
theory  was  the  first  to  include  the  modular  ratio,  n.  Also  about  the 
same  time,  v.  Thullie  (ref.  7)  in  1897  and  Ritter  (ref.  8)  in  1899 
published  the  first  theories  on  the  ultimate  strength  of  reinforced 
concrete  which  considered  the  inelastic  or  nonlinear  stress-strain 
relationship  for  concrete  in  the  compression  block.  Ritter  was  the 
first  to  introduce  the  parabolic  stress  distribution.  However,  most  of 
the  investigators  of  this  period  agreed  that  the  Coignet/Tedesco  straight 
line  method  was  accurate  enough  for  design  purposes  and  its  simplicity 
was  appealing  from  a computational  viewpoint.  Hence,  their  theory 
became  the  ;."i  theory  at  the  turn  of  the  century  and  led  to  rapid 

developments  in  the  use  of  reinforced  concrete  as  a construction  material 

After  the  turn  of  the  century  and  the  general  acceptance  of  the  straight- 
line  method  of  design,  only  a few  investigators  continued  to  research 
ultimate  strength  and  inelastic  stress  distribution.  Assumptions  that 
were  generally  common  to  all  of  the  proposed  new  theories  included  the 
validity  of  Bernoulli's  hypothesis  regarding  strain  proportional  to  the 
distance  from  the  neutral  axis,  no  bond  slip  between  the  concrete  and 
reinforcing  steel,  and  the  concrete  carried  no  tension.  The  reinforcing 
steel  was  assumed  to  have  a trapezoidal  stress-strain  relationship. 
Normally,  two  modes  of  failure  were  considered: 

(1)  Tension  Fai lure--after  the  reinforcing  steel  yields,  the 
concrete  crushes  in  compression  as  a result  of  the  upward 
movement  of  the  neutral  axis. 

(2)  Compression  Failure--the  concrete  crushes  in  compression 
prior  to  yielding  of  the  tensile  reinforcement. 
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V.  ►.(lewn  P.  Neumann  (1890)  ■ ..iqneT  and  -ed.  -..  ; 1 14 

’ '‘'^inddi  j Theory  ' 
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f.  V.  t’  *.erijer  (p.'j-i  f.  ,ufnscn  G.  v.  Kdjir.czy  (1*^33) 

^ . i.  ' ''net'  A.  Rrdn<1tjdei;  (193S) 

' r . S.  rthitnev  (1037) 


I . J.  Mensch  (1914) 


[after  Hognestad  (ref.  2)] 

Figure  i.  Flexural  Analysis  Theories  (1  of  2) 


The  few  investigators  that  continued  with  the  research  on  ultimate 
strength  included  Talbot  (refs.  9 and  10),  Emperger  (ref.  11),  Suenson 
(ref.  12),  Mensch  (ref.  13),  Kempton  Dyson  (ref.  14),  and  Lyse,  Slater,  and 
Zipprodt  (refs.  15  and  16).  Emperger  and  Suenson's  proposed  theories  were 
based  on  a rectangular  stress-block  (fig.  1)  but  considered  tension  failures 
only.  Hence,  no  compatibility  equation  was  needed.  Force  equilibrium 
required  that 


f"ab  = A 
c 


s^ 


(1) 


where 

f“  = flexural  strength  of  concrete  in  compression 
a = depth  of  stress-block 
b = width  of  cross  section 
= area  of  tensile  reinforcement 
f^  = stress  in  tensile  reinforcement 

Moment  equilibrium  implied  that 

<2> 

where 

M = moment  capacity  of  section 

d = distance  from  compressive  face  of  section  to  centroid 
of  reinforcing  steel 


By  combining  eqs. 
where 


(1)  and  (2)  and  substituting  the  reinforcement  ratio,  p. 


the  moment  capacity  can  be  expressed  in  the  form  that  became  well  known  in 
later  years;  i .e. , 


bd^” 


(3) 


Suenson  let  f^  equal  f^^  (cube  strength)  and  f^  equal  the  yield  strength  of 
the  reinforcing  steel. 


Mensch's  proposed  theory  (ref.  13)  in  1914  was  based  on  Ritter's  parabolic 
stress  distribution  and  included  both  tension  and  compression  failures. 
Compression  rei nforcement  could  also  be  included  in  the  design  method.  No 
assumption  was  made  regarding  ultimate  strain  in  the  concrete  nor  was  compat- 
ibility used.  For  tension  failures,  the  computations  were  very  similar  to 
Suenson's;  for  compression  failures,  he  assumed  that  the  neutral  axis  was  lo- 
cated at  the  centroid  of  the  tension  reinforcement.  This  resulted  in  the 
following  moment  capacity: 


where  f^  is  the  compressive  strength  of  a test  cylinder  or  prism.  Mensch 
considered  this  limiting  condition  too  extreme,  however,  and  suggested 
for  t'ei'ifor  "it  >.  ■ a factor  of  1/2.6.  Mensch  also  pointed  out 

that  the  standard  theory  did  not  agree  with  test  results  at  high  loads. 


Kempton  Dyson  (ref.  14)  in  1922  proposed  a theory  based  on  an  elliptical 
stress  distribution.  His  results  were  not  significantly  different  from 
Mensch's.  For  compression  failures  he  found  that 


During  the  first  two  decades  of  this  century,  the  standard  theory  became 
so  widely  used  that  its  approximate  nature  was  often  forgotten.  The 
allowable  concrete  compressive  stress,  f^,  was  normally  taken  as  0.325f^, 
where  f^  was  the  strength  of  a 6-by-12-in  cylinder.  From  this  it  was  normally 
concluded  that  the  safety  factor  against  compression  failure  was  about  three. 
During  this  period,  however,  Slater  and  Zipprodt  (ref.  15)  and  Slater  and 
Lyse  (ref.  16)  pointed  out  that  the  actual  safety  factor  was  much  larger  than 
the  ratio  this  reemphasized  the  actual  inelastic  behavior  of  concrete. 


Prior  to  1920,  bending  stresses  were  ignored  in  the  design  of  reinforced  con- 
crete columns.  It  was  assumed  that  these  stresses  were  provided  for  in  the 
safety  factor.  These  stresses  were  not  considered  mainly  because  suitable 
structural  analysis  methods  for  continuous  structures  were  not  available. 


This  lack  was  eliminated  by  the  development  of  the  Slope  Deflection  Method 
(ref.  17)  and  later  by  the  Moment  Distribution  Method  (ref.  18). 

Another  significant  development  took  place  in  the  early  1920s  when  McMillan 
(ref.  19)  published  a study  on  column  test  data  that  showed  that,  because  of 
creep,  reinforced  concrete  columns  under  a load  may  develop  steel  stresses 
much  higher  than  those  predicted  by  the  straight-line  method.  McMillan's 
study  led  to  an  extensive  ACI  column  investigation  in  the  1930s  by  Lyse, 
Slater,  and  Richart,  which  was  reported  by  Richart  and  Brown  (ref.  20). 

Their  work  resulted  in  the  development  of  equations  for  the  ultimate  strength 
of  axially  loaded  columns. 

Another  milestone  in  the  advancement  of  reinforced  concrete  occurred  in  1931; 
Emperger  (ref.  21)  published  a critical  study  of  the  modular  ratio  and  allow- 
able stresses. 

Prompted  in  the  United  States  by  the  ACI  column  investigation  and  in  Europe 
by  Emperger's  paper,  a renewed  interest  in  the  ultimate  strength  behavior  of 
reinforced  concrete  began  in  the  early  1930s.  In  1932,  Stiissi  (ref.  22) 
published  a paper  in  which  he  discussed  the  ultimate  strength  of  beams  with 
tension  reinforcement  only.  Considered  in  his  study  were  compression  and 
tension  failures  and  the  effect  of  strain  hardening  in  the  reinforcement. 
Brittle  failures  (failures  due  to  rupture  of  the  tension  steel  immediately 
after  the  formation  of  tension  cracks  in  the  concrete)  were  also  discussed 
in  terms  of  small  percentages  of  reinforcing  steel;  Stiissi 's  method  of 
analysis  of  the  ultimate  moment  capacity  was  very  general  and  several 
theories  developed  later  were  actually  refinements  and  improvements  on 
his  work.  Stiissi 's  theory  was  based  on  an  arbitrary  form  of  the  compres-.ion 
stress-block  but  characterized  by  the  stress-block  constants  ki  and  k , the 
compressive  strength  in  flexure,  f^,  and  an  ultimate  strain,  ■ (fig.  1). 

The  constants  kj  and  k?  relate  to  the  magnitude  and  location  of  the  internal 
compressive  force  in  the  concrete.  The  compressive  force  is  therefore  l>;b  t 
He  assumed  f|(  was  equal  to  f^.  Stiissi  determined  his  slress-bloik  constant- 
from  concentric  compression  tests  on  prisms.  He  determined  k-  - 0.70  to  0.77 
k?  = 0.39  to  0.41,  and  = 0.0020  to  0.0025.  These  parameters  were  smaller 
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than  those  reported  by  later  investigators  because  he  did  not  recognize  that 
larger  strains  could  be  developed  in  bending. 

From  the  assumptions  in  figure  1,  equilibrium  of  forces  requires  that 

A f = kibcf 
s s ‘ c 

where  c is  the  distance  from  the  compressive  face  of  the  member  to  the 
neutral  axis. 

Moment  equilibrium  yields 

M = kibcf^(d  - k2c) 

For  tension  failures  with  f^  equal  to  f^  (the  yield  strength  of  the  rein- 
forcement), the  ultimate  moment  capacity  of  a section  was  expressed  as 


Stussi  suggested  that  kr/ki  = 0.55  was  accurate  enough  for  practical  purposes. 
For  compression  failur-es,  using  the  straight-line  strain  distribution,  he  de- 
veloped the  following  compatibility  equation: 


where 

= modulus  of  elasticity  of  steel 
s 

^ = steel  strain 

One  of  Stussi 's  conclusions  was  that  the  safety  factor  with  the  standard 
theory  was  2.3  to  4.1. 


Schreyer  (ref.  23)  in  1933  reported  on  a Stussi -type  theory  for  which  he  es- 
tablished the  stress-strain  relationship  for  concrete  in  compression  from 
cube  tests.  The  stress-strain  curve  was  hyperbol ical ly  shaped  with  ^ 0.0063. 

His  value  for  was  independent  of  concrete  stremth  and  his  resulting  ex- 
pressions for  moment  capacity  were,  algebraical ly,  quite  involved. 
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The  ultimate  moment  capacity  of  rectanqular  beams  failing  in  tension  was  dis- 
cussed by  Kazinczy  (ref.  24)  in  1933.  His  theory  was  based  on  a rectanqular 

stress-block.  It  was  assumed  that  f = f and  f"  = f ' ; therefore,  the  re- 

s y c c 

suiting  expression  of  eq.  (3)  was  used  for  the  moment  capacity  of  the  beans. 

A major  portion  of  Kazinczy's  work  consisted  of  the  analysis  of  test  results 
of  two-span  continuous  beams--one  of  the  first  papers  dealing  with  a plastic 
theory  for  reinforced  concrete  beams. 

In  1934,  Gebauer  (ref.  25)  reported  on  a theory  for  rectangular  beams  failing 
in  tension.  His  theory  was  also  based  on  a rectangular  stress-block;  however, 
he  included  some  tensile  strength  of  the  concrete.  Gebauer  assumed  that  the 
tension  was  the  result  of  shrinkage  of  the  concrete  around  the  reinforcement. 

Also  in  1934,  Baumann  (ref.  26)  reported  on  a study  of  the  buckling  of  rein- 
forced concrete  columns  subjected  to  concentric  or  eccentric  loads.  In  his 
study  a relationship  between  moment  and  rotation  for  reinforced  concrete  mem- 
bers was  required.  Baumann  found  that  a parabola  satisfactorily  approximated 
the  stress-strain  curve  for  concrete.  However,  he  discovered  through  tests  of 
eccentrically  loaded  prisms  that  the  ultimate  strain  in  flexure,  c^,  was  larger 
than  the  corresponding  strain  at  the  maximum  concrete  stress,  c^,  for  a 
concentric  compression  test.  For  a concrete  strength  of  3500  psi,  he  deter- 
mined ^ = 0.0018  and  ^ = 0.0025  to  0.0033.  Baumann,  therefore,  used  the 
stress-block  shown  in  figure  1. 

Bittner  published  two  papers  (refs.  27  and  28)  in  1935  and  1936  on  the  in- 
elastic behavior  of  reinforced  concrete.  For  tension  failures  he  assumed  a 
rectangular  stress-block  and  thus  the  moment  capacity  given  by  eg.  (3)  re- 
sulted. For  compression  failures  he  used  a Stlissi-type  analysis  with  a 
stress-block  similar  to  Baumann's.  However,  Bittner  used  a constant  value 
of  ■ Q = 0.0015,  regardless  of  concrete  strength.  For  his  analysis  he  used 
^ - 0.003,  0.005,  and  0.007  but  did  not  make  a recommendation  as  to  which 
value  should  be  used  in  design. 

Between  1935  and  1937,  Brandtzaeg  (refs.  29  through  32)  reported  on  studies 
that  represented  the  first  complete  analysis  of  the  ultimate  strength  capac- 
ity of  rectangular,  reinforced  concrete  sections.  These  studies  included  the 
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effect  of  compression  rei nforcement  and  bending  in  combination  with  axial 
load.  For  tension  failures  Brandtzaeg  used  a rectangular  stress-block; 
for  compression  failures,  he  used  a stress-block  similar  to  Baumann's.  How- 
ever, it  was  improved  by  the  development  of  a relationship  between  the  ulti- 
mate strain,  and  the  compressive  strength  of  the  concrete.  Also,  f^  was 
assumed  equal  to  0.85f\  Brandtzaeg  also  introduced  a ■ 'r.'?*'*'',  (n  = 

into  his  expressions  for  ultimate  moment.  He  verified  his  theory  by 
comparing  the  results  to  the  results  from  tests  of  20  beams,  13  eccentrical ly 
loaded  columns,  and  several  auxiliary  specimens. 

Emperger,  who  in  1931  had  written  the  paper  critical  of  the  modular  ratio  and 
the  allowable  stress  approach  to  the  design  of  reinforced  concrete  members 
(ref.  21),  published  a report  (ref.  33)  in  1936  based  on  his  review  of  5 
years  of  discussion  on  ultimate  strength  design  and  concluded  that  satis- 
factory results  for  the  ultimate  strength  of  reinforced  concrete  beams  could 
be  obtained  with  the  assumptions  shown  in  figure  1. 

In  1936,  Saliger  (ref.  34)  presented  a thorough  study  of  rectangular  beams  in 
which  he  considered  tension,  compression,  brittle,  and  balanced  modes  of  fail- 
ure. (Balanced  failure  is  simultaneous  crushing  of  the  concrete  and  yielding 
of  the  tensile  reinforcement.)  Saliger  used  the  same  basic  approach  as  Stiissi 
(ref.  22),  but  he  assumed  that  ko  = l/2ki;  this  is  equivalent  to  replacing 
the  curved  stress-block  with  a rectangular  one.  He  also  assumed  f^  = f^; 
this  results  in  eq.  (3)  becoming  the  expression  for  the  ultimate  moment  ca- 
pacity of  a beam  failing  in  tension.  The  value  of  kj  was  then  determined  by 
observing  the  position  of  the  neutral  axis  at  failure;  he  found  values  from 
0.90  to  0.94. 


Whitney,  one  of  the  better-known  researchers  in  reinforced  concrete  in  the 
United  States,  published  papers  on  his  ultimate  strength  theories  from  1937 
to  1948  (refs.  35  through  38).  A rectangular  stress-block  was  used  in  the 
analysis  of  tension  failures  with  the  assumption  that  P = 0.85f^.  The 
following  expression  was  then  developed  for  the  ultimate  moment  capacity: 
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Whitney  made  no  assumptions  regarding  ultimate  strains,  strain  distribution, 
or  bond  slip.  For  compression  failures  he  assumed  from  test  results  a limit- 
ing value  of  a/d  = 0.537.  This  led  to  the  following  expression  for  ultimate 
moment : 


JL 

bd^ 


Whitney  also  developed  expressions  for  the  ultimate  strength  of  rectangular 
and  round  sections  subjected  to  axial  load  and  bending  moment  and  failing  in 
either  tension  or  compression. 


In  1941,  Cox  (ref.  39)  reported  on  tests  of  110  rectangular  beams.  Using  a 
rectangular  stress-block,  he  developed  expressions  for  tension  and  compres- 
sion failures  with  and  without  compression  reinforcement.  He  assumed  f^  = 
f^;  this  results  in  eq.  (3)  for  tension  failures.  For  compression  failures, 
a value  for  the  critical  reinforcement,  p^^,  was  experimentally  developed. 


This  results  in 


M ^ J 

bd.  - 2.76  'c 


The  1/2.76  factor  is  between  the  1/2.6  that  Mensch  (ref.  13)  used  and  the 
1/3  that  Whitney  found. 


Jensen  (refs.  40  and  41)  in  1943  published  a very  complete  report  on  rec- 
tangular beams  with  tension  reinforcement  only.  His  analysis  was  of  the 
Stiissi  type,  but  he  considered  a trapezoidal  stress-block.  He  developed 
the  following  equation  for  the  ultimate  moment  capacity  of  a beam: 


where  N was  a function  of  concrete  strength.  For  tension  failures,  N = 2 was 
assumed;  this  is  the  same  as  eq.  (3).  For  compression  failures,  a compati- 
bility equation  was  used. 
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In  1951,  Hognestad  (ref.  2)  published  results  of  an  extensive  investigation 
of  the  ultimate  strength  behavior  of  reinforced  concrete  members  subjected  to 
combined  bending  and  axial  loads.  He  tested  120  specimens,  which  included 
both  square  tied  and  cylindrical  spiral  columns.  In  the  investigation  he  de- 
veloped an  inelastic  flexure  theory  of  the  Stiissi  type.  Hognestad  used  the 
concrete  stress-strai n relationship  shown  in  figure  1.  He  also  compared  his 
results  to  the  theories  of  Whitney  and  Jensen.  Hognestad's  paper  included  a 
compilation  of  previous  theories  and  his  theory  seemed  to  include  the  best 
features  of  all  of  them. 

In  the  early  1950s,  a consolidation  of  information  concerning  ultimate 
strength  design  was  initiated.  Existing  theories  were  extended  to  include 
combined  bending  and  axial  load  as  well  as  prestressed  concrete  members. 
Ultimate  strength  design  methods  were  also  introduced  into  the  building  codes 
of  several  countries.  In  October  1955  an  ASCE-ACI  Joint  Coiiinittee  on 
■ • published  its  final  report  (ref.  42),  which  culminated  more 

than  10  years  of  committee  work.  This  report  led  to  changes  in  the  ACI  Build- 
ing Code,  318-56  (ref.  43),  and  thus  permitted,  for  the  first  time  in  the 
United  States,  the  use  of  ultimate  strength  design  methods  for  reinforced 
concrete  flexural  members. 

After  the  publication  of  the  ASCE-ACI  Joint  Committee  report  and  the  inclusion 
of  ultimate  strength  provisions  in  the  1956  ACI  Building  Code,  Hognestad  (ref. 
44)  published  the  results  of  a review  of  current  literature  regarding  the  in- 
elastic stress  distribution  in  flexure  of  reinforced  concrete  members.  These 
results  were  compared  with  the  recommendations  of  the  Joint  ASCE-ACI  Committee, 
and  from  the  comparison  it  was  concluded  that  the  committee's  design  coeffi- 
cients were  well  substantiated  by  test  results  and  that  a simplified  rectangu- 
lar stress  distribution  gave  satisfactory  accuracy  for  common,  practical  de- 
sign cases. 

In  1961,  Mattock,  Kriz,  and  Hognestad  (ref.  45)  presented  the  development  of 
an  ultimate  strength  design  theory  based  on  an  equivalent  rectangular  stress 
distribution  in  the  concrete  compression  zone.  The  proposed  theory  was  in 
general  accord  with  the  appendix  to  the  1956  ACI  Building  Code,  but  it  had  a 


much  broader  application.  The  method  was  applied  to  a wide  variety  of  rein- 
forced concrete  beams  and  columns  under  various  combinations  of  axial  load 
and  flexure.  These  results  were  compared  to  results  of  experimental  investi- 
gations and  the  agreement  was  excellent.  Since  a wide  range  of  variables 
was  considered  in  this  investigation,  it  was  concluded  that  the  theory  pre- 
dicted ultimate  strength  with  sufficient  accuracy  for  all  types  of  structural 
sections  encountered  in  structural  design,  including  odd-shaped  sections. 

This  work  was  the  basis  for  the  ultimate  strength  provisions  for  flexure  and 
axial  load  in  the  1963  ACI  Building  Code  (ref.  46). 

With  the  general  acceptance  of  ultimate  strength  methods  for  the  design  of 
reinforced  concrete  structures  and  the  acceptance  of  inelastic  behavior  of 
reinforced  concrete,  increased  interest  and  attention  were  directed  toward 
limit  design  of  reinforced  concrete  structures.  Even  at  the  present  time, 
in  standard  design  procedures,  an  inelastic  stress  distribution  is  consid- 
ered in  designing  cross  sections  to  resist  moments  and  loads  determined 
from  an  elastic  analysis. 

In  November  1964,  an  International  Symposium  on  .■T  •nr)  ;-,- 

;■  / '-njyete  (ref.  47)  was  held  to  present  recent  work  directed  speci- 

fically toward  the  goal  of  a more  basic  understanding  of  the  flexural  behav- 
ior of  reinforced  concrete  and  the  elimination  of  the  basic  contradiction  in 
design  philosophy,  e.g.,  ultimate  strength  design  from  an  elastic  analysis. 
Most  of  the  papers  presented  were  concerned  with  some  aspect  of  limit  design 
(e.g..  Mattock's  paper  (ref.  48), 

/It- ' navete  Bciiru' , and  Roy  and  Sozen's  paper  (ref.  49),  rujtilitu  v' 

' K't'c:,:.)  Most  of  the  papers  were  concerned  with  the  behavior  of  reinforced 
concrete  beams  and  slabs  in  the  plastic  region  up  to  maximum  load.  However, 
Barnard  presented  one  of  the  first  reports  (ref.  50)  on  beam  behavior  beyond 
maximum  load  to  collapse.  Based  on  the  concept  of  concrete  as  a 
• ■ material,  Barnard  showed  that  a beam  could  continue  to  rotate  when 

the  bending  moment  was  falling  off,  but  it  would  not  collapse  unless  an  en- 
ergy balance  in  the  beam  ceased  to  be  satisfied. 

A comprehensive  annotated  bibliography  (ref.  51)  which  covers  limit  design 
investigations  between  1917  and  1968  was  prepared  by  Cohn. 
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The  most  recent  investigation  into  the  behavior  of  reinforced  concrete  beams 
in  the  range  beyond  maximum  moment  capacity  was  reported  by  Iqbal  and  Hatcher 
(ref.  52)  in  1975.  The  behavior  of  reinforced  concrete  beams  in  the  post- 
crushing  range  (where  moment  is  decreasing  while  deformations  continue  to 
increase)  was  studied.  A theoretical  model  of  the  failure  mechanism  in  the 
post-crushing  region  was  presented.  The  results  were  compared  with  the  ex- 
perimental results  of  tests  on  six  beams  without  compression  or  transverse 
rei nforcement. 
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SECTion  III 


EXPERIMENTAL  PROGRAM 


The  experimental  phase  of  this  investigation  consisted  of  statically  testing 
to  collapse  17  rectangular,  reinforced  concrete  beams.  The  duration  of  the 
beam  tests  varied  from  2 to  12  minutes.  All  beams  had  the  same  span  length, 
cross  section,  and  hinged  end  supports.  The  two  test  parameters  considered 
were  the  axial-to-lateral-load  ratio  and  the  shear-span- to-beam-depth  ratio. 
For  a two-point  symmetrical  lateral  load,  the  shear  span  is  the  distance  from 
the  support  to  the  load  point.  Figure  2 shows  the  general  loading  scheme  for 
the  beams.  The  axial-to-lateral-load  ratio,  P/F,  remained  constant  for  the 
duration  of  each  test.  The  shear-span  parameter  is  expressed  as  the  a^/d 
ratio;  a^  is  the  shear  span  and  d is  the  distance  from  the  compression  face 
of  the  concrete  to  the  centroid  of  the  tensile  reinforcement  (effective 
depth).  The  parameters  considered  in  this  i nvestiqation  were  a^/d  ratios  of 
5,  4,  and  3 and  P/F  ratios  of  3,  2,  and  0.  Each  beam  test  was  identified  by 
a three-number  designation  (e.g.,  5-3-1).  The  first  number  refers  to  the 
nominal  a^/d  ratio,  the  second  refers  to  the  nominal  P/F  ratio,  and  the  third 
to  the  beam  number  of  that  particular  configuration.  Table  1 presents  the 
test  designations  for  each  configuration. 


TEST  APPARATUS 

The  loading  device  used  in  this  investigation  was  a test  frame  used  by  Crist 
(ref.  53)  in  a reinforced  concrete  deep  beam  study;  but  it  was  modified  to 
accommodate  the  axial  load  application.  Figure  3 shows  the  modified  test 
frame  which  consisted  of  an  upper  portion  that  provided  reaction  for  the  lat- 
eral load  and  a lower  portion  that  provided  the  axial  load  and  support  sys- 
tem. The  two  portions  were  tied  together  by  five  vertical  structural  T- 
sections. 

Lateral  load  was  applied  by  a 100,000-lb-capaci ty  hydraulic  ram  with  a 13- 
in  stroke.  Axial  load  was  applied  by  two  200,000-lb-capacity  double  acting 
rams  with  26-in  strokes  mounted  in  a horizontal  position.  A single  hydraulic 
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Ficjure  3.  Test  Frame 


system  was  used  to  activate  the  rams  to  insure  a constant  P/r  ratio  through- 
out each  test.  The  hydraulic  system  had  a capacity  of  9,500  psi  and  could 
be  operated  from  a portable  remote  control  with  an  adjustable  load  rate. 

The  total  lateral  load  was  divided  into  a two-point  load  by  a steel  distri- 
bution beam  which  imparted  force  to  the  beams  through  2-1/2-in-diameter  rol- 
lers and  4-by-9-by- 3/4-in  steel  bearing  plates.  One  end  of  the  distribution 
beam  was  free  to  translate  and  rotate  while  the  other  end  was  only  free  to 
rotate.  The  bearing  plates  were  seated  to  the  beams  with  a thin  layer  of 
high-strength  gypsum  compound.  The  length  of  the  shear  span  was  established 
by  the  position  of  the  lateral  loads. 

Axial  load  was  applied  to  the  beams  by  using  the  two  horizontal  rams  in  ten- 
sion; this  resulted  in  compression  on  the  beam  because  of  the  pivoting  of  the 
vertical  reaction  arms  (fig.  4).  The  ratio  of  axial  to  lateral  load  was  de- 
fined by  the  pull-point  position  of  the  rams  on  the  reaction  arms.  To  adjust 
the  axial-to-lateral-load  ratio,  the  connection  point  of  the  tension  rams  to 
the  pivot  arm  was  changed.  For  zero  axial  load,  the  horizontal  rams  were  me- 
chanically and  hydraulically  disconnected  from  the  system.  The  hinged  condi- 
tion at  the  ends  of  the  beams  was  insured  by  transmitting  the  axial  load 
through  self-aligning,  roller-bearing  pillow  blocks  and  4-in-diameter  steel 
shafts.  The  axial  force  was  applied  through  the  ■•Za;-:,' ' • acntv  i of  the  beam 
cross  section.  The  plastic  centroid  of  a section,  as  defined  by  the  1963  ACI 
Building  Code  (ref.  46),  is  the  centroid  of  resistance  to  load  computed  under 
the  assumption  that  the  concrete  is  uniformly  stressed  to  f|(  and  the  rein- 
forcing steel  is  uniformly  stressed  to  f^. 

To  insure  that  the  forces  measured  in  the  horizontal  rams  during  the  tests 
could  be  accurately  converted  to  axial  forces  in  the  test  beams,  a series  of 
calibration  tests  was  conducted  prior  to  actual  beam  testing.  The  calibra- 
tion tests  consisted  of  loading  a dummy  beam  axially  while  measuring  both 
the  forces  in  the  hydraulic  rams  with  force  links  and  the  axial  load  in  the 
beam  with  a load  cell.  Under  the  assumption  that  all  joints  in  the  mechani- 
cal linkage  between  the  horizontal  rams  and  the  dummy  beam  were  frictionless, 
the  calculated  axial  load  in  the  beam  was  compared  to  the  load  indicated  by 
the  load  cell.  The  agreement  between  calculated  and  measured  forces  in  the 
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dummy  beam  was  within  10  percent;  however,  the  corrections  derived  from  the 
calibration  tests  were  used  to  determine  axial  loads  in  the  beams  during  the 
actual  tests. 


BEAM  SPECIMENS 
Geometry 

The  beam  span  was  12  ft,  6-1/4  in  measured  from  center  to  center  of  the 
pivotal  shafts  at  the  beam  ends.  This  length  included  the  steel  end  reaction 
devices  shown  in  figures  4 and  5.  Cross-sectional  properties  and  geometry 
were  the  same  for  all  specimens.  The  beams  were  15  in  in  overall  depth  and 
9 in  wide  with  a depth  from  the  compressive  face  of  the  concrete  to  the  cen- 
troid of  the  tensile  steel  of  12-1/2  in.  Tensile  reinforcement  consisted  of 
three  No.  6 (3/4-in-diameter)  bars.  Although  it  was  intended  that  the  speci- 
mens be  essentially  singly  reinforced,  two  No.  2 (1/4-in-diameter)  bars 
were  placed  in  the  top  of  the  beams.  These  were  included  to  assist  in  beam 
fabrication  and  to  facilitate  making  strain  measurements  in  the  compression 
zone  of  the  beams.  No.  2 stirrups  were  placed  at  6-in  intervals  along  the 
length  of  the  beams  (fig.  5). 

To  facilitate  axial  load  application,  the  concrete  portion  of  the  beams  was 
terminated  at  end  bearing  plates  to  which  the  end  reaction  devices  were 
bolted.  The  longitudinal  reinforcement  was  welded  to  the  end  bearing  plates 
to  assure  adequate  anchorage  for  the  bars  and  to  assure  development  of  the 
full  flexural  and  shear  capacities  of  the  beams.  Additional  reinforcement 
was  also  welded  to  the  end  plates  to  provide  a mechanism  for  shear  transfer 
between  the  concrete  and  the  end  supports.  Figure  5 shows  the  details  of 
the  end  bearing  plates. 

Reinforcing  Steel 

The  principal  longitudinal  reinforcing,  which  consisted  of  three  No.  6 
bars,  had  a yield  strength  of  62,000  psi  and  conformed  to  ASTM  Specification 
A615-60.  All  the  steel  was  produced  from  the  same  heat  to  insure  consistency 
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among  the  beam  specimens.  Tensile  tests  were  performed  on  samples  from 
several  lengths  of  the  rebar  to  insure  control  and  to  determine  the  mechan- 
ical properties  of  the  steel.  A typical  stress-strain  curve  is  shown  in 
figure  6a. 

The  stirrups  and  compression  reinforcement  were  intermediate-grade  steel 
conforming  to  ASTM  Specification  A15  with  a yield  strength  of  52,000  psi . 
Although  not  covered  by  Specification  A305,  beams  with  the  fio.  2 bars  had 
deformations  similar  to  those  with  the  No.  6 bars.  A typical  stress-strain 
curve  is  shown  in  figure  6b. 

Concrete 

The  concrete  used  in  the  beams  had  a nominal  compressive  strength  of  5,000 
psi  and  was  produced  with  Type  III  Portland  cement  and  a maximum-size 
aggregate  of  3/8  in.  The  coarse  aggregate  was  mostly  well-rounded  natural 
material  of  uniform  gradation  (100  percent  by  weight  passing  the  3/8-in 
sieve  with  less  than  12  percent  by  weight  passing  the  No.  4 sieve)  with 
less  than  10  percent  by  weight  crushed  material.  The  fine  aggregate  was 
a washed  material  conforming  to  ASTM  Specification  C-33  for  concrete  fine 
aggregate  and  had  less  than  2 percent  by  weight  passing  the  No.  200  sieve. 

The  fineness  modulus  of  the  sand  varied  from  2.6  to  3.1. 

The  concrete  was  mixed  at  the  Eric  H.  Wang  Civil  Engineering  Research  Facility 
(CERE)  in  a 16-ft  nontilting,  rotating  drum,  electric-powered  concrete  mixer. 
Three  to  six  control  cylinders  for  each  beam  were  cast  in  6-in-diameter,  12- 
in-high,  waxed  cardboard  molds.  The  beams  were  cast  in  their  normal  position 
in  steel  forms.  Each  beam  and  its  control  cylinders  were  cast  from  a single 
batch  of  concrete.  The  beam  concrete  was  compacted  with  a 1-in-diameter,  12- 
in-long,  electric  vibrator  probe  which  operated  at  about  10,000  rpm.  The  con- 
trol cylinders  were  compacted  by  vertical  vibration  at  a frequency  of  10  Hz 
for  about  1-1/2  minutes  on  a vibrating  table  designed  and  fabricated  at  CERE. 

The  beams  and  control  cylinders  were  cured  under  polyethylene  plastic  sheets 
for  at  least  48  hours  before  the  forms  were  removed  and  the  cylinder  molds 
stripped.  The  beams  and  cylinders  were  then  left  to  cure  together  under  the 
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polyethylene  sheets  until  about  2 weeks  prior  to  the  beam  test,  at  which 
time  curing  continued  at  ambient  conditions. 


The  control  cylinders  were  tested  to  failure  the  day  of  the  beam  test.  Ten- 
sile splitting  tests  and  stress-strain  curves  were  made  for  a limited  number 
of  the  control  cylinders.  Figure  7 shows  a typical  stress-strain  curve  for 
the  concrete.  The  portion  of  the  curve  beyond  maximum  stress  could  not  be 
measured  because  of  the  characteristics  of  the  testing  machine  used.  Table  2 
presents  the  concrete  strengths  for  the  various  beams. 

INSTRUMENTATION 

Instrumentation  for  all  beams,  with  a few  exceptions,  was  similar.  Measure- 
ments were  made  at  nodal  locations  that  corresponded  to  nodes  used  in  the  be- 
havioral model  development.  Figure  8 shows  the  location  of  the  reference 
nodes.  Measurement  locations  were  symmetrical  about  the  beam  centerline. 

Measurements  made  included  vertical  deflections  along  the  beam,  vertical  and 
horizontal  deflections  at  the  support  shafts,  relative  horizontal  displace- 
ments along  the  beam,  steel  and  concrete  strain  at  various  locations,  rota- 
tions at  the  beam  ends,  and  lateral  and  horizontal  loads. 

Displacement 

Vertical  deflection  measurements  were  made  at  seven  locations  along  the  beam 
(1,  2,  4,  and  6 in  figure  8).  At  the  beam  centerline,  station  6,  the  meas- 
urement was  made  from  the  ground  to  the  bottom  surface  of  the  beam.  At  the 
other  six  locations,  the  measurements  were  made  from  the  ground  to  aluminum 
brackets  bonded  to  the  beam  side  at  middepth.  Support  movement  during  the 
tests  was  monitored  by  vertical  and  horizontal  displacement  measurements  on 
both  sides  of  the  support  shafts  at  each  end  of  the  beams.  1 inear  potenti- 
ometers with  gage  lengths  from  1 to  6 in,  depending  on  the  location  of  the 
measurement , were  used  for  these  measurements. 
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Table  2.  Test-Day  Compressive  Strength  of  Concrete 


Standard  deviations  are  shown  in  parenthesis 


Horizontal  displacement  measurements  were  also  made  at  six  locations  along 
the  beams.  These  measurements  were  made  as  relative  displacements  between 
stations  1 and  2,  2 and  4,  and  4 and  6.  The  measurements  were  made  with  di- 
rect current  differential  transformer  (DCDT)  displacement  transducers  with 
full-scale  ranges  of  ♦ 0.05  in. 

Strain 

Steel  strain  in  the  longitudinal  tensile  reinforcement  was  measured  at  sta- 
tions 2 through  6.  The  measurements  were  made  on  the  middle  reinforcing 
bar.  Strain  measurements  on  the  two  compression  steel  bars  were  also  taken 
at  stations  2 through  6,  except  under  the  load  points,  in  which  case  no  meas- 
urement was  made.  At  station  6 (the  beam  centerline)  a strain  gage  was 
mounted  on  each  bar.  At  the  remaining  stations,  the  measurements  were  taken 
on  alternate  sides  of  the  beams.  Strain  measurements  were  also  taken  in  the 
first  two  stirrups  outside  the  load  points  in  all  but  the  first  four  beams 
tested,  in  which  the  gages  were  inadvertantly  left  out.  Measurements  of  the 
steel  strain  were  made  with  350-ohm,  epoxy-backed,  foil  strain  gages  which 
had  a 1/2-in  gage  length  and  a gage  factor  of  2.125. 

The  steel  strain  gages  were  mounted  with  epoxy  cement  on  a widened  and 
smoothed  portion  of  a longitudinal  rib  of  the  bars.  Lead  wires  were  then 
attached  and  the  gage  was  waterproofed  with  a plastic  sealant.  The  lead 
wires  exited  the  beams  through  holes  in  the  side  forms. 

All  concrete  strain  measurements  were  made  on  the  beam  surface.  Three  con- 
crete strain  measurements  were  taken  at  the  beam  centerl ine--one  at  the 
middle  of  the  top  surface,  one  1 in  from  the  top  on  the  side  of  the  beam, 
and  one  2-1/2  in  from  the  top  on  the  side.  In  addition,  in  the  5-series 
(i.e.,  a^/d  = 5),  the  two  side  measurements  were  also  made  at  station  4; 
in  the  4-series,  the  two  side  measurements  were  also  made  at  stations  5 
and  3;  and  in  the  3-series  the  side  measurements  were  also  made  at  stations 
5 and  2. 

Concrete  strain  measurements  were  made  with  300-ohm,  paper-backed , wire 
strain  gages  which  had  a 1-in  gage  length  and  a gage  factor  of  2.05.  The 
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surface  on  which  a strain  gage  was  to  be  mounted  was  first  ground  or  sanded 
to  remove  foreign  material  and  to  smooth  the  surface.  The  surface  was 
cleaned,  filled  with  epoxy,  and  sanded  smooth.  Epoxy  was  then  used  to  bond 
the  gage  to  the  surface. 

Rotation 

Rotation  measurements  were  made  at  both  ends  of  the  beams  with  gages  spe- 
cially fabricated  at  CERF.  The  rotation  gages  consisted  of  a pendulum  sus- 
pended from  the  paddle  portion  of  a OX  type  velocity  gage.  The  rotation  of 
the  pendulum  relative  to  the  gage  body  was  measured  by  the  variable  induc- 
tance transducer  of  the  velocity  gage. 

Data  Acquisition  and  Reduction 

The  electrical  instrumentation  measurements  were  continuously  recorded  on 
1-in  magnetic  tape.  Recording  was  accomplished  at  a tape  speed  of  30  in/ 
sec;  thus,  each  test  could  have  a maximum  recording  time  of  15  minutes. 
However,  to  reduce  the  computer  time  associated  with  the  digitizing  of  the 
analog  tapes,  it  was  intended  that  each  test  be  less  than  7-1/2  minutes. 
Analog  data  tapes  were  digitized  and  data  reduction  and  presentation  were 
performed  at  the  Kirtland  Air  Force  Base  computing  facilities. 

Photoelastic  Coating 

In  addition  to  the  electrical  measurements,  photoelastic  coating  was  used 
on  the  sides  of  the  beams  to  determine  the  complete  concrete  strain  distri- 
bution in  compression  in  the  constant  moment  region.  Since  only  the  com- 
pression strain  pattern  was  of  interest,  the  coating  was  applied  to  the 
top  half  of  the  beams  only. 

The  strain  pattern  was  observed  by  taping  a Polaroid  sheet  over  the  coating 
and  then  photographing  the  resulting  fringe  patterns  that  occurred  during 
loading  of  the  specimens.  Slow-speed  (six  frames  per  second),  color  iiwtion 
pictures  were  taken  of  the  coated  area  during  the  tests.  Theoretically, 
based  on  the  characteristics  and  thickness  of  the  coating,  the  number  and 
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order  of  the  fringes  can  be  related  to  the  strain  in  the  coating.  To  relate 
the  st»ain  pattern  to  a load  condition,  a digital  voltneter,  which  displayed 
the  output  from  the  lateral  load  force  link,  was  located  on  the  test  frame 
and  photographed  with  the  strain  patterns. 


SECTION  IV 

DEVELOPMENT  OF  ANALYTICAL  MODEL 


A behavioral  model  for  static  response  prediction  of  hinge-ended  reinforced 
concrete  beams  under  combined  axial  and  lateral  loads  was  formulated.  The 
beam  response  was  analyzed  by  dividing  the  beam  along  its  length  into  a num- 
ber of  segments  connected  at  points  called  :e.-  and  assuming  the  curvature 
varied  linearly  between  these  nodes.  The  lateral  loads  were  applied  as  point 
loads  at  the  nodes  and  the  axial  loads  were  applied  longitudinally  at  the  end 
nodes  (hinge  points).  Only  beams  of  symmetrical  geometry  and  loads  about  the 
centerline  were  considered.  A constant  axial-to-lateral-load  ratio  was  main- 
tained throughout  the  calculated  bea;:i  response.  A maximum  concrete  strain 
was  not  specified  in  the  model;  consequently,  the  collapse  point  was  not  de- 
termined. However,  the  response  was  described  beyond  any  significant  load- 
carrying  capacity.  A CDC  6600  digital  computer  was  used  to  solve  for  the 
beam  behavior  (i.e.,  deflections,  concrete  and  rebar  strains,  and  curvatures 
at  the  node  points  as  a function  of  applied  loading)  from  the  developed  model. 

Results  of  the  behavioral  model  calculation,  plotted  with  the  experimental 
results,  are  presented  in  appendix  A.  Appendix  B presents  the  computer  pro- 
gram used  for  the  analytical  calculations. 

MATERIAL  BEHAVIOR 

The  behavioral  model  used  predicted  the  response  well  beyond  the  maximum 
resistance  of  the  beam.  It  was  necessary,  therefore,  that  the  concrete 
stress-stra'i n relationship  used  in  this  model  development  represent  the  com- 
plete stress-strain  behavior  and  not  just  that  to  maximum  stress.  When  a 
standard  universal  testing  machine  is  used  to  test  concrete  cylinders,  the 
specimens  normally  crush  when  the  maximum  stress  is  reached.  This  behavior 
leads  to  difficulties  in  obtaining  data  on  the  stress-strain  relationship 
of  concrete  beyond  the  point  of  maximum  stress.  Barnard  (ref.  50)  and  other 
researchers  cited  by  Barnard  have  tested  concrete  cylinders  using  specially 
constructed  very  stiff  constant-strain-rate  testing  machines.  These  machines 


permitted  concrete  specimens  to  be  tested  to  very  large  compressive  strains, 
well  beyond  those  occurring  at  maximum  stress.  Figure  9 shows  the  complete 
stress-strai n curve  for  concrete  obtained  with  the  special  testing  machines. 
Based  on  the  descending  portion  of  the  curve  where  stress  decreases  with  in- 
creasing strain,  Barnard  describes  concrete  as  a .ti-  r ■ material 

(unlike  mild  steel  which  is  a r-  i-i!  material).  This  strain- 

softening characteristic  of  concrete  has  a significant  effect  on  the  behav- 
ior of  the  beam  after  maximum  resistance  has  been  attained. 

Figure  10  shows  the  concrete  stress-strain  curve  used  for  this  model  formu- 
lation. This  curve  was  assumed  to  represent  the  complete  stress-strain  be- 
havior of  the  concrete  in  the  beam.  The  initial  portion  of  the  curve  is  de- 
fined by  the  parabola 


where 

f = concrete  stress 
c 

f^  = flexural  strength  of  concrete  in  compression 
= concrete  strain 

c = strain  at  maximum  concrete  stress 
0 

This  parabolic  form  was  first  used  by  Hognestad  (ref.  2)  and  suosequently 
by  many  other  investigators. 


The  second  portion  of  the  curve  is  a descending  straight  line  connecting 

the  top  of  the  parabola  at  a strain  of  ^ and  a horizontal  line  at  a 

stress  of  0.2f". 

c 

The  third  portion  of  the  stress-strain  curve,  a horizontal  line,  suggests 
that  concrete  can  sustain  a stress  of  0-2f^  to  infinity.  This  approach 
has  been  used  previously  by  Barnard  (ref.  54),  Yamashiro  and  Siess  (ref. 
55),  and  Kent  and  Park  (ref.  56). 


The  maximum  moment  capacity  of  a beam  with  no  axial  load,  singly  reinforced 
with  steel  having  a bilinear  stress-strain  relationship,  and  the  concrete 


strain,  in/in  x 10"’ 

[after  Barnard  (ref.  50)] 

Figure  9.  Typical  Complete  Stress-Strain  Curve  for  Concrete 
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Figure  10.  Concrete  Stress-Strain  Curve  Used  in  Analytical  Model 


stress-strain  characteristics  shown  in  figure  10  occurs  at  a strain  between 
■ 0 and  • 0 (the  strain  at  which  the  concrete  stress  becomes  20  percent  of 
f|l).  The  exact  value  of  the  strain  corresponding  to  the  maximum  moment  ca- 
pacity depends  on  the  slope  of  the  descending  straight  line  and  occurs  at 
the  strain  where  the  ratio  of  the  area  under  the  stress-strain  curve  to  the 
area  of  the  rectangle  f^  , designated  ki,  becomes  a maximum.  The  strain  at 
which  ki  becomes  maximum  is  designated  as  • Figure  11  further  illustrates 
the  relationship  between  the  concrete  stress-strain  curve  and  kj. 


With  reference  to  figure  11  and  for 


the  area.  A,  under  the  concrete  stress-strain  curve  is 


A = f" 
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The  ratio  of  the  area  under  the  concrete  stress-strain  curve  to  the  area 
of  the  rectangle  f^  is  k;. 
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The  strain  at  which  ki  is  maximum  can  be  determined  from 
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By  designating  ^ the  strain  corresponding  to  maximum  ki,  the  slope  of  the 
descending  portion  of  the  curve  required  to  make  k maximum  at  is 


Why  the  maximum  moment  capacity  of  an  under-reinforced  concrete  flexural 
member  (tension  failure)  occurs  at  maximum  kj  can  be  explained  by  considering 
the  internal  forces  in  the  beam.  After  the  tension  reinforcement  yields 


and  the  curvature  is  increasing,  the  tension  force,  T^,  remains  constant. 

Force  equilibrium  requires  that  the  compression  force  in  the  concrete, 

C^,  be  equal  to  the  tension  force,  or 

C = T 
c s 

However , 

C,  = k f"  be 
C 1 c 

where  b is  the  beam  width  and  c is  the  distance  from  the  compression  face 

of  the  beam  to  the  neutral  axis.  Therefore,  with  increasing  curvature  prior 

to  a concrete  strain  of  , k increases;  this  requires  a decrease  in  c. 

0 1 

As  the  neutral  axis  moves  upward,  the  internal  moment  arm  between  and 
also  increases  and  this  results  in  an  increasing  moment.  Beyond  an  ex- 
treme fiber  strain  of  ^ and  as  the  curvature  increases,  k decreases. 
Consequently,  c must  increase;  this  results  in  a decreasing  internal  moment 
arm  and  decreasing  moment. 
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The  tensile  strength  of  the  concrete  was  not  ignored.  Concrete  cracking  was 
specified  by  a limiting  concrete  tensile  strain,  The  modulus  of  elas- 

ticity of  concrete  in  tension  was  assumed  as  suggested  in  the  ACl  Building 
Code  ( ref.  46) ; i . e. , 


E 


cone 


i-n  psi 


where  w is  the  unit  weight  of  concrete  (145  Ib/ft').  Unloading  of  the  con- 
crete was  assumed  to  occur  at  a slope  parallel  to  the  initial  slope  of  the 
compression  portion  of  the  curve. 
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Figure  12  illustrates  the  stress-strain  relationship  assumed  for  the  tension 
reinforcing  steel.  The  curve  is  bilinear  , nd  includes  a second-degree  curve 
representation  of  strain  hardening.  The  equation  of  the  curve  beyond  strain 
hardening  for  >.  _ is 


f 

s 


\h 


^sh^  ■ '^sh^^  ■ 


(4) 


where  f^  is  the  steel  stress  and  f^  is 


the  yield  stress. 
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Figure  12.  Typical  Stress-Strain  Curve  for  Reinforcing  Steel 


■ = steel  strain 

, = strain  at  start  of  strain  hardening 
sh 

A ^ = initial  slope  of  strain  hardening  curve 

B , = constant  determined  to  match  actual  steel  stress-strain  behavior 
sh 

In  compression,  the  rei nforcement  v/as  assumed  to  behave  as  in  tension  except 
with  no  strain  hardening.  Unloading  of  the  steel  was  assumed  to  occur  at 
the  initial  slope  of  the  curve,  E^. 

MODEL  FORMULATION 

A coimion  method  for  solving  the  response  of  a statically  loaded  beam  is  to 
increment  the  applied  load  and  calculate  the  resulting  strains,  stresses,  ro- 
tations, and  deflections.  This  method,  however,  presents  numerical  instabil- 
ity problems  when  the  material  involved  has  a strain-softening  characteristic 
and  the  response  is  calculated  into  the  decreasing-load  region.  Conseguently , 


the  method  used  to  calculate  the  responses  in  this  investigation  consisted 
of  incrementing  the  strain  at  the  top  of  the  beam  at  the  centerline  and 
then  computing  the  associated  loads,  rotations,  deflections,  and  the  strains 
and  resulting  stresses  at  the  centerline  and  at  the  remaining  node  points. 
The  solution  was  repeated  until  the  strain  at  the  top  of  the  beam  reached 
some  large  value  (e.g.,  0.05). 


BEAM  BEHAVIOR 

Figure  13  illustrates  the  model  used  for  the  overall  response  calculations. 
The  axial  load,  P,  and  the  total  lateral  load,  F,  on  half  the  beam  are  re- 
lated by  the  factor  K;  i.e., 

K = P/F 

For  no  axial  load  (i.e.,  K = 0),  a very  small  K was  used  in  the  calcula- 
tions (e.g.,  0.005)  because  subseguent  calculations  involved  division  by  K. 


Figure  13.  Beam  Model  for  Overall  Response 
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Distribution  of  the  lateral  load  was  specified  by  the  factors  A.,  i = 1,  N, 
such  that 

f.  = A.F  = A.P/K,  0 _ A.  _ 1 .0  (i 

The  shear  force  in  any  segment  and  the  bending  moment  at  any  node  were  ex- 
pressed in  terms  of  the  axial  load,  P,  by 


V.  = V.P,  i=l,  N 


m.  = M.P,  i=l , N 


where 


V.  = shear  force  in  segment  i 

nm  = bending  moment  at  node  i due  to  lateral  load  only 


Vi  = 1.0/K 


^ = ^--1  - 


Mo  = 0.0 


M.  = . V.C.,  i=l,  N 


Node  deflections  and  rotations  can  be  calculated  if  curvatures,  are 
known  at  each  node  and  the  distribution  of  the  curvature  between  nodes 
is  assumed  to  be  linear.  The  curvatures,  were  determined  from  force 
eguilibrium  at  the  nodes.  It  was  also  assumed,  then,  for  the  calculation 
of  the  deflections,  that  the  beam  rotations  were  concentrated  at  the  nodes 
and  segments  between  the  nodes  remained  straight.  The  angle  change  or 
change  in  rotation  at  a node  i was  computed  as  follows: 

AO.  = (2^..  + /6 

The  slope  of  the  segment  at  the  hinged  end  was  computed  as  follows: 
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^ 

n 

= ^ 

M. 

i = l 

>< 

The  slope,  then,  of  any  segment  was 

i-1 

0.  = 

0i  - 

■,e . 

1 

j=i  ' 

The  beam  deflection  at  any  node  was 

computed  as  follows: 

^•-1 

+ O.f. 

i 

1 1 

wi  th 

0 

0 

BEAM  SECTION  EQUILIBRIUM 

Figure  14  illustrates  the  stresses  and  associated  forces  acting  on  a beam 
section  resulting  from  a known  strain  distribution.  The  assumptions  asso- 
ciated with  section  equilibrium  are  as  follows; 

(1)  Strain  is  proportional  to  the  distance  from  the 
neutral  axis. 

(2)  Concrete  and  steel  stress-strain  relationships 
are  as  described  previously. 

(3)  Shear  behavior  does  not  affect  the  flexural 
behavior  of  the  beams. 

The  force  resultants  shown  in  figure  14  are  discussed  below. 


Compressive  Concrete  Force  (C^) 


C^  was  determined  by  integrating  the  concrete  stress-strain 
cur've  presented  in  figure  10. 


k.  - ^ L-i 
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bkdk  1 


where 

b 

kd 

k 


k 


C 


c 


bkdk I k 


beam  width 

distance  from  compressive  face  of  beam  to  neutral  axis 
factor  relating  flexural  strength  of  concrete  in  beam 
to  concrete  cylinder  strength 
strain  in  concrete  at  extreme  fiber  of  bean 
concrete  stress-strain  curve  shape  factor  previously 
di scussed 


Force  in  Compression  Reinforcetnent  (C^) 


a;  f; 


where 

A'  = area  of  compression  reinforcement 
s 

f^  " stress  in  compression  reinforcement 

However,  when  ' , the  compression  rei nforccment  is  assumed 

s - u 

to  buckle  and  no  longer  carry  load. 

Tensile  Force  in  Concrete  (T^) 

T = i F bX 

c 2 cone  cr  c 

where 

= depth  of  tension  stress  block  (determined  from 
strain  distribution) 


Tensile  Force  in  Main  Reinforcement  (T^) 


T 

s 


A 


s 


f 

s 
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where 


= area  of  reinforcinq  steel 
= stress  in  reinforcing  steel 


Horizontal  force  equilibrium  requires  that 


P = C ^ C - T 


The  axial  force,  P,  acts  through  the  of  the 

beam  section  which  is  located  a distance  Xp^  from  the  bottom 
of  the  beam.  The  plastic  centroid  of  a section  is  the  centroid 
of  resistance  to  load  computed  under  the  assumption  that  the 
concrete  is  uniformly  stressed  to  k f^  and  the  reinforcing  steel 
is  uniformly  stressed  to  f . 

y 


Moment  equilibrium  about  the  plastic  centroid,  Xp^,  yields  an  expression 
for  the  moment  resistance  of  the  beam. 


Mp  = C^(moment  arm  abt  PC)  C^(moment  arm  abt  PC) 
t T^(moment  arm  abt  PC)  + T^(moment  arm  abt  PC) 


RESPOriSC  SOLUTION 

The  concrete  strain  at  the  top  of  the  centerline  section  was  incremented 
and  the  associated  forces  and  displacements  were  calculated.  The  follow- 
ing steps  were  employed  in  determining  the  response: 

(1)  A set  of  displacements  (• .,  i=0,  M)  was  assumed  at  the 

start  of  each  response  cycle.  For  the  first  cycle,  the 
deflections  were  assumed  to  be  zero  i = 0,  N)  and 

subsequently  the  final  deflections  calculated  from  the 
previous  cycle  were  used  as  the  assumed  deflections  for 
the  next  cycle. 

(2)  With  at  node  N)  and  the  load  ratio  (K  = P/F) 

known,  the  strain  at  the  bottom  of  the  section  was  incre- 
mented until 


49 


where 


m 


+ P\ 


= moment  resistance  of  the  section  at  node  N 
[eq.  (7)] 

m^j  = moment  at  node  h due  to  applied  lateral  loads 
[eq.  (6)] 

= assumed  deflection  at  node  N 


At  this  time  in  the  solution,  an  axial  load,  P,  and  the 
corresponding  moments  at  each  node  were  known  as  well  as 
the  curvature  at  the  centerline,  :|,j. 

(3)  With  the  axial  load  and  moments  known  at  the  remaining 
nodes,  the  concrete  strain  and  the  strain  at  the  bottom 
of  the  beam  at  each  remaining  node  were  incremented 
until  the  resisting  and  applied  forces  were  equal;  i.e., 

Mj,.  = m.  + PA.  i = l,  ri-1 

K 1 ^ 1 


This  established  a set  of  curvatures  ( : , i = l , W-1). 

(4)  With  the  new  set  of  curvatures  known,  a revised  sot  of 
nodal  deflections  was  calculated.  The  new  deflection 
at  the  centerline  was  compared  to  the  previously  assumed 
deflection.  If  the  agreement  was  within  a specified 
tolerance,  the  cycle  was  ended  and  a new  cycle  was  started 
with  the  next  increment  of  at  node  N.  If  the  center- 
line  deflections  did  not  agree  with  the  specified  tolerance, 
the  old  deflections  were  replaced  with  the  new  and  the  cycle 
was  repeated  with  the  same  ^ at  node  .N.  The  tolerance  used 
in  comparing  new  and  previous  deflections  was  a relative 
error  of  0. 01 , or 


'S  - A , . 
new  old 


/A 


new 


0.01 


There  was  no  failure  (collapse)  criterion  associated  with  the  model;  there 
fore,  the  response  was  calculated  to  some  large  ^ (e.a.,  0.05).  The  term 
inal  value  of  ^ was  selected  large  enough  to  obtain  the  response  well 
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beyond  the  maximum  beam  resistance  and  far  enough  to  establish  significant 
behavioral  characteristics. 

Beyond  the  maximum  moment  resistance  of  a beam,  the  curvature  at  the  node 
of  maximum  moment  (the  centerline  node)  continues  to  increase  vnth  decreas- 
ing moment.  However,  at  the  remaining  nodes  the  curvature  and  associated 
strains  decrease  with  decreasing  moment.  As  the  beam  is  forced  through 
further  deflection,  the  curvature  at  the  center  becomes  more  concentrated 
as  the  remaining  portions  of  the  beam  unload  and  decrease  in  curvature. 

To  adequately  model  the  unloading  of  portions  of  the  beam,  the  maximum 
strains  encountered  at  each  node  were  stored  for  comparison  with  subsequent 
strains.  Also,  based  on  observed  behavior  beyond  maximum  moment,  the  por- 
tion of  beam  undergoing  the  extremely  large  curvature  has  a finite  length 
that  must  be  accounted  for  in  the  deflection  calculations.  Therefore,  at 
the  center  node,  whenever  the  strain  at  the  compression  face  was  greater 
than  c^,  deflections  were  calculated  based  on  an  assumed  width  of  maximum 
curvature  of  d (the  beam  depth). 


ANALYTICAL  MODEL  PARAMETERS 

The  geometric  parameters  used  were  the  same  as  those  for  the  experimental 
beam  specimens  described  in  section  3 and  are  summarized  in  table  3. 

Table  4 presents  the  compressive  concrete  strength,  f^,  and  the  load  ratio, 

K = P/E,  used  in  the  calculations.  The  concrete  strengths  presented  are 

the  average  of  several  (3  to  6)  compression  tests  performed  the  day  of  the 

beam  tests  on  6-by-12-in  cylinders  cast  with  the  beams.  The  load  ratios 

were  determined  from  lateral  and  horizontal  for^e  measurements  made  during 

«►  • 

the  beam  tests.  ", 

The  load  factors,  A^ , that  relate  the  individual  concentrated  nodal  forces 
to  the  total  lateral  load  on  half  of  the  beam  (eg.  5)  were  as  follows: 
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Table  3.  Geometric  Parameters 


Parameter 

Descri ption 

Value 

Beam  Width 

9.0  in 

d 

Depth  of  Tension  Rei nforcement 
from  Compressive  Face  of  Beam 

12.5  in 

d' 

■ 

Depth  of  Compressive  Reinforce- 
ment from  Compressive  Face  of 
Beam 

1 . 5 i n j 

1 

j 

t 

Total  Beam  Depth 

15.0  in 

A 

Area  of  Tensile  Reinforcement 

L ^ 

1.32  in 

i 

1 

-H 

j 

A' 

s 

Area  of  Compressive  Reinforce- 
ment 

0.10  in 

L 

Beam  Length 

150.25  in 

N 

Number  of  Nodes 

. 

1 

th 

Length  of  i Beam  Segment 

1 1 

= 10.625  in 

1 

= 16.5  in 

i 

L_ 

i' , 

. ..  r = 12.0 

s ^ • 

Table  4.  Average  Concrete  Compressive 
Strength  and  Load  Ratio 


Beam 

Designation 

f , 
c 

- Jill  . 

Load  Ratio  ^ 

(K  = P/F) 

5-0-1 

5014 

0.005 

5-0-2 

4980 

0.005  1 

5-2-1 

4336 

1.91  1 

5-2-2 

4651 

1.88 

5-3-2 

4753 

3.10 

5-3-3 

5583 

3.21 

4-0-1 

5288 

0.005 

4-2-1 

5205 

1.97 

4-2-2 

4990 

1.88 

4-3-1 

5205 

3.18 

4-3-2 

5423 

3.15 

3-0-1 

5142 

0.005 

3-2-1 

4792 

1.85 

3-2-2 

4831 

1.87 

3-3-1 

5252 

3.03 

3-3-2 

5028 

3.08 

3-3-3 

4344 

3.11 

Scries  5 


^ 1 . < u f.  ~ 0 

^ 1 » ^ 

As  = 1.0 

Series  4 

A 1 3 = 0 

> ' ^ * 

Au  = 1.0 

Series  3 

A)  , . s ( - 0 

^ » 

A<  = 1.0 

Table  5 presents  the  parameters  associated  with  the  concrete  model.  The 
modulus  of  elasticity,  used  only  when  the  tensile  strenqth  of  the  concrete 
was  considered,  was  the  value  suggested  in  the  ACl  Buildirif)  Code.  The  un- 
load modulus  was  arbitrarily  selected  as  the  initial  tanqent  to  the  para- 
bolic portion  of  the  stress-strain  curve.  The  strain  at  maximum  stress, 
was  selected  as  an  average  from  the  cylinder  tests  on  the  beam  concrete. 


Table  5.  Parameters  Associated  With  Concrete  Model 


Parameter 

Descri ption 

Value 

Modulus  of  Elasticity 

57,400  'Y?^  psi 

C 1 , 

unload 

Modulus  of  Elasticity  for 
Unloading  Concrete 

2f 

_c 

0 

0 

Strain  at  Maximum  Stress 

0.0022 

1 ^ 

Strain  at  Maxinum  k 

1 

0.0035 

cr 

strain  at  Concrete  Cracking 

0.0001  i 

k 

1 

1 

Ratio  of  Flexure  Strength 
of  Concrete  to  Cylinder 
Strength 

1.00 

The  few  tensile  splitting  tests  conducted  on  the  concrete  indicated  an  ap- 
proximate tensile  strength  of  400  psi.  This  value  with  a modulus  of  4 x 10' 
psi  leads  to  the  value  of  the  cracking  strain,  that  was  used  in  the 

calculations.  The  value  of  l used  was  selected  based  on  some  observed  be- 

u 

havior  and  also  on  a value  that  seemed  to  yield  better  comparisons  in  some 
areas  between  calculated  and  experimental  beam  response. 

The  tensile  reinforcement  was  assumed  to  be  elastic-plastic  and  to  include 
strain  hardening.  The  strain  hardening  portion  of  the  steel  behavior  was 
described  by  a second-degree  curve  [eg.  (4)]. 

The  compression  reinforcement  was  also  assumed  to  be  elastic-plastic,  but 
not  to  include  strain  hardening.  Also,  the  compression  reinforcement  was 
assumed  to  have  no  strength  beyond  a strain  of  t:^.  Table  6 presents  the 
reinforcing  steel  parameters  that  were  determined  from  tensile  tests  on 
the  rebar. 

The  calculations  of  the  beam  behavior  were  made  to  a maximum  compressive 
strain  at  the  top  of  the  beam  of  0.05  in/ in. 


Table  6.  Reinforcing  Steel  Parameters 


! Parameter 

r-  - 

E 

U- 

Description 

Value 

Modulus  of  Elasticity 

30  X 10'  psi 

r 

f 

y 

Yield  Stress  for  Tensile 
Reinforcement 

' 

— 

62.0  ksi 

f ' 

y 

Yield  Stress  for  Compressive 
Reinforcement 

52.0  ksi 

sh 

Strain  at  Onset  of  Strain 
Hardening  (Tensile  reinforce- 
ment only) 

0.003 

\h 

Parameters  Used  to  Describe 
Stra  in-Hardening  Portion  of 
Tensile  Rei nforcemen L 

2."^  X 10’  psi  ' 

1 

2^.6  X 10'  psi  ' 

J 

SECTION  V 
RESULTS 


Results  from  both  the  experimental  investigation  and  the  behavioral  model 
analysis  are  presented.  Analytical  data  were  obtained  for  the  node  points 
at  which  experimental  measurements  were  made.  When  possible,  both  results 
are  presented  together  so  that  comparisons  may  be  made  and  the  results  dis- 
cussed. Appendix  A presents  all  tiie  measured  experimental  and  the  calcu- 
lated analytical  data  for  each  beam. 


GENERAL  BEHAVIOR 

The  general  response  of  the  beams  can  be  illustrated  by  their  load-centerl i ne 
deflection  curves.  Fiqurr  15  presents  an  idealized  response  curve.  The  be- 
havior of  the  beams  was  generally  stiff  for  the  initial  10  percent  of  the 


f iqur'e  15.  Tyoical  l.oad-Centerl  i ne  Deflection  Curve 
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loadinq,  followed  by  a fairly  linear  portion  of  reduced  stiffness.  Upon 
yieldinq  of  the  tensile  reinforcement  the  curve  flattened  considerably  un- 
til the  'iiaximum  load  was  achieved.  The  behavior  beyond  yield  was  a func- 
tion of  the  maqnitude  of  the  axial  load.  When  the  axial  load  was  zero,  or 
very  small,  the  flattened  portion  of  the  curve  was  relatively  lon().  With 
hiqh  axial  load  the  flattened  portion  was  either  short  or  nonexistent. 

Upon  reachinq  maximum  load,  with  no  axial  load,  the  load  decreased  at  a 
qradual  rate  with  increasinq  deflection.  With  hiqh  axial  load,  the  load 
dropped  rapidly  upon  reaching  the  maximum  load.  Figure  16  presents  the 
load-centerline  deflection  curves  for  the  17  beams  tested. 


MODE  OF  FAILURE  AMD  CRACK  PATTERN 

The  mode  of  failure  of  the  beams  was  flexural  tension;  i.e.,  the  tension 
reinforcement  began  yielding  before  the  concrete  crushed.  The  beams  col- 
lapsed when  the  concrete  at  the  top  of  the  beams  was  crushed.  Based  on 
the  ultimate  shear  strength  of  the  beams  predicted  by  the  ACI  Building 
Code,  the  beams  with  short  shear  spans  could  have  failed  in  shear.  How- 
ever, the  small  percentage  of  shear  rei nforcement  was  very  effective  iti 
preventing  shear  failure.  In  the  analytical  calculations,  shear  behavior 
was  not  considered;  therefore,  only  flexural  failures  were  predicted. 

Also,  prediction  of  crack  spacing  and  width  was  beyond  the  scope  of  the 
analytical  effort;  consequently,  only  observed  crack  behavior  is  reported. 

Because  the  experimental  test  measurements  were  continuously  recorded  on 
magnetic  tape,  the  duration  of  each  test  could  be  no  longer  than  15  min- 
utes (maximum  time  on  one  roll  of  tape  when  recording  at  30  in/sec).  How- 
ever, to  conserve  magnetic  tape  and  to  reduce  data-reduction  costs,  the 
test  times  were  held  to  2 to  12  minutes.  These  relatively  short  test 
times  precluded  the  marking  of  cracks  at  specific  load  increments,  which 
is  commonly  done  in  static  tests  of  reinforced  concrete  components. 

Following  each  test,  the  cracks  were  marked  and  the  beams  were  photographed. 
Figure  17  presents  the  final  crack  patterns.  Visible  crack  formation  was 
documented  by  slow-speed  motion  pictures  taken  to  record  the  fringe  patterns 
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from  the  photoelastic  coating.  The  motion  pictures  indicated  that  crack- 
irui  first  occurred  in  the  middle  of  the  constant-moment  region  of  the 
beams  and  then  spread  tov-jard  the  supports.  The  initial  cracks  generally 
formed  at  the  stirrups.  At  advanced  stages  of  loading,  additional  cracks 
formed  between  the  initial  cracks.  Near  the  maximum  load  capacity  of  the 
beams,  horizontal  cracks  formed  at  the  level  of  the  tensile  reinforcing  steel; 
this  indicated  bond  failure.  In  general,  the  beams  with  high  axial  load 
were  cracked  less  than  those  with  no  axial  load.  It  can  also  be  seen  from 
ti<)ure  1/’  that  inclined  diagonal  tension  cracks  extended  from  the  tensile 
re i tiforcement  level  to  the  vicinity  of  the  load  points.  Table  7 presents 
the  percentages  of  maximum  load  at  which  first  visible  cracking,  first 
diauonal  cracking,  and  concrete  spalling  occurred  (based  on  the  slow-speed 
motion  pictures). 


Table  7.  Concrete  Cracking  and  Spalling  Data 


Load,  Percentage  of  Maximum 


Beam 

Designation 

t First 

1 Flexural 
1 Cracking 

First  Visible 
Diagonal  Cracking 

Concrete 

Spalling 

5-0-1 

' 14 

1 39 

95 

5-2-2 

i 12 

34 

100 

5-3-3 

23 

38 

100 

4-0-1 

12 

32 

92 

1 

4-2-2 

i 16 

64 

100 

4-3-1 

25 

74 

100 

3-0-1 

15 

1 

1 

; 52 

100 

3-2-1 

13 

1 59 

100 

3-3-2 

22 

; ..__8i 

100 

Note:  These  values  were  taken  from  motion  pictures 

of  the  tests. 
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FLEXURAL  BEHAVIOR 


The  agreement  between  calculated  and  measured  behavior  prior  to  cracking 
of  the  concrete  was  good.  The  agreement  from  cracking  to  maximum  load 
resistance  of  the  beams  is  suiiaiiarized  in  table  8.  Comparisons  are  made 
for  yield  and  maximum  loads  and  for  deflections  at  yield  and  maximum 
loads.  Load  agreement  was  good.  Comparisons  ranged  from  -3  to  9 percent 
for  yield  load  and  from  -7  to  9 percent  for  maximum  loads,  with  an  aver- 
age in  both  cases  of  3 percent.  The  deflection  agreement  was  not  as  good 
as  the  load  agreement  (generally  the  case  in  reinforced  concrete  investi- 
gations). All  measured  yield  deflections  were  higher  than  those  calcu- 
lated; this  indicated  that  the  beams  were  not  as  stiff  as  predicted. 

Yield  deflection  ratios  varied  from  10  to  36  percent,  with  an  average  of 
24  percent.  Deflection  agreement  at  maximum  load  was  much  more  erratic 
than  the  yield  deflection  agreement;  even  though  the  average  was  -2 
percent,  the  ratios  ranged  from  -21  to  44  percent. 

Deflections  at  collapse  or  the  actual  collapse  load  could  not  be  compared 
since  the  theoretical  calculations  did  not  include  a criterion  for  this 
phenomenon.  The  calculations  were  arbitrarily  continued  to  a compression 
strain  at  the  top  of  the  beam  of  0.05  in/in.  The  actual  collapse  behavior 
of  a beam  is  related  to  the  stiffness  of  the  testing  apparatus  and  the  en- 
ergy released  into  the  crushing  portion  of  the  beam.  Therefore,  collapse 
behavior  in  this  experimenta  configuration  would  not  be  applicable  to  pro- 
totype behavior  under  real  loads. 

A comparison  of  the  tensile  steel  strain  data  for  the  various  beams  (ap- 
pendix A)  shows  that  the  agreement  between  calculated  and  measured  strain 
became  worse  with  distance  from  the  load  point.  This  indicates  that  bond 
failure  or  slip  between  the  end  of  the  beam  and  the  load  point  is  at  least 
partially  the  reason  why  the  actual  deflections  prior  to  yield  were  greater 
than  the  calculated  deflections. 
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Table  8.  Summary  of  Calculated  arid  Measured  Load  and  Centerline  Deflection  Data 


Standard  Deviation  |0. 03310. 055  0.072 


DUCTILITY 


A common  way  to  express  the  ductility  of  reinforced  concrete  beams  is  to 

compare  deflection  of  the  beam  at  first  yielding  of  the  tensile  steel, 

Ayieid>  to  the  deflection  at  maximum  load,  Table  9 presents  both 

calculated  and  measured  ductility  ratios  (;i  = A... /A  . , ,)  and  measured 

collapse  ductility  ratios  = A ,,  /A  . , where  A ,,  is  the 
^ ' collapse  yield'  collapse 

deflection  of  the  beam  at  collapse.  The  ductility  ratios  decreased  as 
the  axial  load  increased.  This  occurs  because  the  higher  axial  stress  in 
the  compression  zone  of  the  beam  causes  failure  of  the  concrete  prior  to  a 
significant  amount  of  tensile  steel  yielding  and  accompanying  deflection. 


HINGE  FORMATION 

Normally,  hinge  formation  is  associated  with  beam  behavior  after  yielding 
of  the  tensile  reinforcement  and  before  crushing  of  the  concrete.  Beams 
subjected  to  a single  concentrated  load  or  fixed-end  beams  display  very 
pronounced  hinge  formations  at  tiie  concentrated  load  or  at  the  fixed  ends. 
In  these  cases  the  hinge  is  concentrated  at  the  maximum  moment  point  and 
is  confined  to  a finite  length  of  the  beam  because  of  the  moment  qradient. 
In  this  investigation,  because  of  the  two-point  loading,  which  results  in 
a constant-moment  region,  this  type  of  hinging  was  not  experienced,  even 
though  there  was  a small  moment  gradient  in  some  of  the  beams  because  of 
the  P-  effect.  There  was  instead  a general  yielding  of  the  tensile  rein- 
forcement along  the  constant-moment  region.  However,  a secondary  hinge 
formed  in  the  constant-moment  region.  This  hinge  was  formed  when  the  con- 
crete, at  some  point,  became  stressed  beyond  its  maximum  load-carrying 
capacity  and  entered  the  strain-softening  region.  At  this  point,  the 
load-carry ing  capacity  of  the  beam  began  to  decrease.  Thereafter,  the 
strain  at  the  secondary  hinge  point  continued  to  increase,  while  the 
strain  in  the  remainder  of  the  beam  decreased.  The  farther  the  beam  de- 
flected, the  more  concentrated  the  hinge  became.  Final  collapse  of  the 
beams  occurred  when  the  secondary  hinge  region  could  accommodate  no  more 
rotation  and  almost  completely  disintegrated.  Usually  a diagonal  crack 
formed  across  the  hinge  region  simultaneously  with  the  collapse  of  the 
beam. 
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INTERACTION 


Interaction  between  axial  load  and  bending  moment  was  demonstrated  by  a 
comparison  of  beam  behavior  as  the  axial  load  was  increased  and  also  by 
an  ■’nteraction  diagram  of  axial  load  versus  applied  moment  for  both  cal- 
culated and  measured  behavior.  Figure  18  presents  both  measured  and 
calculated  centerline  deflection  versus  total  lateral  load  for  each  series 
of  beam  tests.  The  increase  in  maximum  load  and  the  decrease  in  ductility 
with  increased  axial  load  are  evident.  Also,  the  variation  in  behavior 
with  a variation  in  concrete  strength  is  demonstrated.  Table  10  presents 
a sui'i'iary  of  interaction  data  at  maximum  load  for  all  beams.  These  data 
are  plotted  on  the  theoretical  interaction  diagram  (fig.  19).  which  was 
derived  with  the  analytical  model  developed  in  section  4. 


ANALYTICAL  MODEL 


In  evaluating  the  results  of  this  investigation,  some  areas  of  possible 
error  in  the  calculated  beam  behavior  should  be  considered.  The  method 
used  to  calculate  the  beam  behavior  involved  replacing  a continuous 
structure  with  an  assemblage  of  discrete  elements  connected  at  nodes. 

Thus,  the  resulting  calculated  beam  response  is  an  approximation  of  the 
actual  behavior.  The  accuracy  of  this  approximation  depends  on  the  ele- 
ment size  and  the  mathematical  representation  of  the  material  character- 
istics of  the  beam.  Also,  in  calculating  the  nodal  rotations  and  deflec- 
tions, the  distribution  of  curvature  between  nodes  was  assumed  to  be 
linear.  Because  of  the  nonlinear  moment-curvature  relationship,  a higher- 
order  assumption  would  result  in  less  error. 


Calculated 


Measured 


Figure  18.  Calculated  and  Measured  Load-Centerline  Deflection  Curves 
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SECTION  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 


CONCLUSIONS 

One  of  the  main  objectives  of  this  investigation  was  to  study  the  static 
behavior  of  reinforced  concrete  beams  under  combined  flexural,  axial,  and 
shear  effects.  However,  none  of  the  beams  failed  in  shear  and,  therefore, 
no  real  data  on  shear  interaction  or  transition  of  failure  modes  under  the 
combined  effects  were  obtained.  The  investigation  did  provide  good  infor- 
mation on  ductility  and  behavior  of  beams  in  the  large  deflection  region. 

In  general,  the  procedures  and  methods  employed  in  this  investigation  for 
beam  specimen  fabrication,  testing,  and  experimental  data  acquisition  were 
carefully  controlled  in  order  to  produce  the  best  results  possible.  Each 
size  of  reinforcement  was  controled  by  using  bars  from  the  same  heat.  The 
concrete,  however,  varied  in  compressive  strength  between  -13  and  12  per- 
cent from  the  nominal  5000  psi.  Although  load  application  techniques  and 
control  were  kept  as  constant  as  possible  for  all  tests,  the  total  time  of 
each  test  varied  from  2 to  12  minutes  (2  minutes  is  not  short  enough  to 
cause  dynamic  effects  and  12  minutes  is  not  long  enough  to  involve  creep). 
But,  a more  constant  test  time  would  have  been  more  desirable. 

Most  of  the  instrumentation  and  the  data-acquisi tion  techniques  used  pro- 
duced satisfactory  results.  There  were  some  exceptions,  however.  The 
data  from  the  strain  gages  mounted  on  stirrups  that  were  expected  to  be 
in  the  region  of  diagonal  tension  shear  cracks  were  very  erratic  and  of 
very  little  value  in  determining  beam  response.  Another  measurement  that 
provided  very  little  insight  into  beam  behavior  was  the  horizontal  meas- 
urement made  between  stations.  The  most  notable  data-acquisi tion  method 
that  did  not  provide  satisfactory  results  was  the  photoelastic  coating. 

This  coating  was  applied  to  obtain  data  on  hinge  formation  in  the  constant- 
moment  region  of  the  beams.  However,  the  nonhomogeneity  of  the  concrete 
caused  irregular  and  erratic  fringe  patterns  and,  consequently,  no  useful 
data  were  obtained.  Also,  near  maximum  load,  the  concrete  began  to  spall. 
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causing  the  coating  to  debond.  The  one  positive  aspect  of  using  the 
photoelastic  material  was  the  slow-speed  motion  pictures  that  provided 
good  visual  documentation  of  the  sequence  of  events  in  the  response  of 
the  beams. 

The  general  beam  behavior  calculated  from  the  analytical  model  agreed 
well  with  the  measured  results,  especially  in  the  region  up  to  maximum 
load.  The  measured  beam  stiffness  (the  slope  of  the  load-deflection 
curve  between  first  concrete  cracking  and  yielding  of  the  tensile  rein- 
forcement) was  always  less  than  that  predicted;  this  seemed  to  be  the 
result  of  bond  slip  or  failure  between  the  rebar  and  the  concrete.  (This 
was  not  considered  in  the  analytical  model.)  The  yield  and  maximum  loads 
predicted  by  the  model  were  within  9 percent  of  the  test  results.  Deflec- 
tion calculations  were  not  as  good  as  load  calculations.  Calculated  yield 
deflecv.ons  were  all  lower  (10  to  36  percent)  than  those  predicted.  Com- 
parison of  calculated  and  measured  deflections  at  maximum  load  showed 
very  erratic  differences.  The  ratio  of  measured  to  calculated  deflections 
ranged  from  -21  to  44  percent. 

The  analytical  model  did  not  calculate  the  collapse  point  of  the  beams; 
it  calculated  the  behavior  to  an  arbitrary  compression  strain  (0.05  in/ 
in)  at  the  top  of  the  beam.  Therefore,  no  comparison  could  be  made  of 
collapse  behavior;  however,  beyond  maximum  load,  the  general  behavior 
predicted  by  the  model  was  good. 


RECOMMENDATIONS 

Since  the  objective  of  this  investigation  was  to  study  beam  behavior  to  beam 
collapse  under  combined  flexural,  axial,  and  shear  forces  and  since  the  mode 
of  failure  for  all  test  beams  was  flexural  tension,  it  is  recommended  that 
another  group  of  beams  without  stirrups  be  tested  under  the  same  general 
loading  conditions.  This  would  provide  information  on  the  transition  from 
a flexural  tension  mode  of  failure  to  a shear  mode. 
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It  is  further  recommended  that  another  group  of  beams  be  tested  with  axial- 
to-lateral-Joaa  ratios  high  enough  to  provide  data  above  the  balance  point 
on  the  axial  load-bending  moment  interaction  diagram.  This  group  should  in- 
clude beams  with  and  without  stirrups  in  order  to  provide  shear- interaction 
da  ta . 


Further  work  with  the  behavioral  model  is  suggested.  There  are  three  ob- 
vious areas  not  considered  in  the  present  model:  (1)  bond  behavior  be- 

tween the  reinforcing  steel  and  the  concrete,  (2)  calculations  of  beam 
shear  strength,  and  (3)  calculation  of  tfie  point  of  collapse  of  the  beam. 
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APPENDIX  A 

CALCULATED  AND  EXPERIMENTAL  DATA 


This  appendix  presents  both  experimental  and  ana- 
lytical results  of  the  behavior  of  the  17  beams 
tested.  Where  possible,  the  corresponding  analyt- 
ical and  experimental  data  are  presented  on  one 
graph  for  comparison. 
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rot. 


NATIONAL  BUREAU  OF  STANDARDS 

MICROCOPY  RESOLUTION  TEST  CHART 


Tal  : jTerjl  Li  .iJ 


TotJl  Lateral  L'>ad.  *ip^  fotdl  Lateral  LooJ. 


Toul  Lateral  Loil,  kips  Total  Lateral  Load,  kips  Total  Lateral  Load,  kips  Total  Lateral  Load,  kips 


proi  rii.JiP-i  1» 


Tor  i|  laTfr4l  L.'fci 


Total  Lateral  load 


APPENDIX  B 


PROGRAM  BEAM 


A listing  of  the  computer  code 
used  to  calculate  the  beam  be- 
havior from  the  analytical 
model  is  presented. 
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O O O 


PIVO'-RA'-'  llEA.  ■ ( !,'.P  JT  . Jo  T PUT  . I A PL  6 . T A Pt  ? . I APLb  , T A^o  y . r I ) 

r 

C ^RO",p.'  . K>-A.'  C-J'iPUlL,'  TiIl  ’JliimvIOP  J-  A S j'.  L 1 I C A o L t Lof^LLo 

C Kt  [ ''tP  ORC  1.  o Co'.L  t.  T i JCA-'  >.  ,'i '(3  I \t.  ly  o'-if'.RAL  A.,ij  AxiAl.j  LohJj. 

JP!;A!LL'  Lit  Cn'-’LtK  ir’.lV?':' 

T L 6 = L A i_  o o'  P ♦ r ^ L ' 1 1 L l\o  f i L ■ • T A P L n = ovjo L lJ 

C-''t‘or;/A/A^t  - *lC'  ^P*f  P»AotTr  t 

i i io“'R  * * r.'  « f f * L Co  ')  * A f,  n,  8 » Af  1 

C *LJl('0)  )to0'  »P*lI  1^0^  ♦oL'J.4lLoJ  »r  ♦ 

: f 1 ( .? 0 1 • -1 . r 0 ) « Lf-  ( ? -\) o > I <:  ~ ) , . ( 2 c ; . ll  1 1 1 2 j ) • t-  1 2 1 2 c < 

C/‘'''JN/:<  / f ”•'  X ( 2^  ' •'Ll'  •■  < ( ) . (,<^7'’Ax  ( 2"|  1 

r ( / L / /,<  1 ( ? c ) ’ * V A 1 2 0 ) 

(■■■'-'. '.Oiw'G/  Al  ( 2C  I 22  I ( 20  I • A.i  ' I 2C  )•  PLL  i A ( 2;)  )•  f ML  i /<  i i o ; »o7ft-  I ^ I <;j 

1 ' . i..-:  ( 20  I » oA  ( 20  ’ • ..  ( 2 ^ - L-P  I ( ) »oLu  I d ( o(i  ) 

f C ■'■,  ( CClI.  , A ) - • ill  ^ o ■>.'•:  I ( ■ . 0*C'o..L,/  ( f ij  i * aAn.o  i , 

1 F0P'-'/u(l-] 

2 ^ OK'-'A  T ( H - IC.  ■'  I 

; CL^  1 

Pr.A'^  ! \L  oT  DA  T A 


C 

c 


2 00 

f-0  0 
‘^''1 

80] 


1 OO'I 

3 

3 

6 

I 

8 


i\  = L>]  (oF^r:  Pv'ii.r.s 

P t A r\  1 , rj 

AL(  I ) = i_>  \oTii  >Jt  r.  L AAi  rP  li-.L!  i\ 

Pt  A^  2 . I AL I 1 1 . I - 1 . N ) 
i r ''  A j I -V  I A L P P Op  r ’ 1 L S 

t^L*'  ' i^tYp 

i<L-»'  2yFtv»L.. ».l^>fl»LbriP«^*LHP*.)J^^P 

Lt  - Or  TK  I CF  L PpoPERi  ILb 
PI  > 2 , I .0.1  .PD  , -lu,  Pv 

I-.  lad  I I A ( .;  ■ L T r ' S 

hLA.  D 2.tAf(l).I  = i.i*) 

iij  AOO  1-].*. 

[I  ( D ( I ) ♦ 1 1 • • ) '“  1 b T ■ o 0 A 

''Or.T  I 

I Apr  I 

I\lAP  IjTA\i_  Lr.A.-'  LL'.uTM.  L^At,  < A j i ij  , C'Jl.CKi  T“  3 T Pl.'.u  Tf1 
Rt  ad  2 . ; L .Ax , XXX .FDf 

i;ilTIAL'2c'  OATA 

DO  1000  1-1*6 
.LlTHI 1 )=0.C 
6A( I ) =0.0 
t C ( I ) =0*0 
f M r ) =0.0 
► FP ( I 1=0. 

THI  T a ( I ) =0.0 

PHI ( I ) =0.0 

FORMAT ( I A . 1 It  12.  ' ) 

FuPPf,  1( // / 36X  X Aa  IAll  V LUADLD  wlA-’  oL  HAV  1 < / / / i 

FOR'-'A  I ( / / 38X  ^bPA'.  A-jL  LoAL  t A|.' A -.l  I LKo  V / ) 

F OR,>' A I ( / / 2 b X ♦ bt  X i I .J  i ' KOF'L  I'  I 1 L j ' / / ) 

F JKMA  r ( / /-3bX  •■  • A I L < I At  P.vv^PtD  I 1 t / / I 


tl 
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^ f T ( / / ' *3  X * 1 *'\/  I ul^  MArv  A r L i\  i S T I 5 **  / / ) 

1?  ' l.  .'„  . A I ( ' ,_l'.  Ml  AA  r ; o-Aa  i Ai-  ■ / 3X  ■ . I \(..- Ic ) / L A ; ,.aAl  - / ^ 3a*uoAu* 

i / ) 

: 3 y JK.V  A T i 1 1 ^ \ P.  ■ ) 
i A f-  ATilA»ltX»  2i  I ^ ^ i 

1>  ATt//-*  '3JuL/‘'llX--'t-0'’L*<T«A^''fr,ACIIvyA  w(»/-  oL\J>’;_i''ii*‘yA'*(-L‘*CjfM*t> 

2x  « MALf  -SPA,\s  / 1 7a  - i I .-.CH'  S ) 3..X  " i-oA'i"  / ) 

16  r jiA'' A T ( * r J I lALt'  I n ; I-  r L ■I.  I i VA  - A ••  A I . . T t-  C.„‘-  t al--j  i o.< » A a ” A i c^c,  ' 'J  a ^ 

Ia-'P  ofdhL  ijHt.MA  jl_*‘'/bA‘(ri'*'‘,  A-;yL^  It*  it  t*OA"i*-)I  — it-  UL^‘  ] I ^'~u  H 

i f FCF  •;  T ACil  a t.A ^ I.  f. t At  .L  AlI.  v l a r a.  ,Lit  / I 

1 7 I t,'.\f,A  I ( 6a  Jt-C  A »u  ..*  1 i A ■ - i Oa  - L Aa  * . ■,  I-  I c ‘ A A*-y  i 1-  - -.^U  ' OA  “ I t ‘ A A ^ f I r A 

1 1 ":x»-'  ''yttaX«LA!fS  VX't  AShA  / / ) 

^ f-  ■■,  T A / X A L -J'AL'  A :>X  - 1,  L F Lt  ^ 1 i ^t,  ■•  I j a * ^ 1 \ A i '.  *•  u X a ^ i r*  A i 1 . * 7 A -'■  ^t  u.  a • o a * _.ml 
A A A * *'  X ♦ * 7 H r "*  R "t 

'•  2 0 X » ‘ I F I ( ui ! tt  / d 1 A ' I V . X C t ■’  T X C I A V,  ! I t v A v i TV-/' 

Id  A ..p'-Aj  ( Itil  ) 


PRIM  r 

r-.l’OT  FATA 

PRIM 

19 

pr:  m 

3 

PRI  \T 

5 

Pi*  ! N T 

12 

PRINT 

1 ' . L . A R 

PAINT 

lb 

print 

1 A . ( ; . A p ( 1 I , ■ n ( i ) , 1 - 1 . , I 

PRIM 

; 

print 

16 

PP  I N T 

l''.  tP.£3.l)t.  .*■  »PV 

PRINT 

A 

P R I ■'I  T 

1 ' 

P i<  I N T 

1 . *A  < _■  * 3 • t,  . r,  t.  R • P r . L X*  * *'’  • t Y * ' . t ^ k ♦ c,  1 1 . A j 1 1 

print 

i 

PRINT 

IS 

COVRJTF  cc,-tCTA-,:r- 

ff''C0‘i  = ‘j5  Af!  . ■ ■ ;=<  T ( t .■•;  t ; aoo  • ) 

AX6=  T /: 

A A 8 = "I  p / r 

roT-3./S;jRT  (FPC*  lAnr.  i 
CP  = P1>C>D  »L 
C ■ 1 = c p » ") 

Ab=FY/FPC-AA  3 
AC  = P YD/F  PC-Ax  3 

XPC  = ( ( Aa  -•  » A*P5«  ' •!  . ) / ^ .■»»-■  • Ak. « t mMj- An3  1 ♦ t • Ao  » ( MAP  - i . I ) / l Aao  *Ak  J + Pr'» 
1 A C + P » A R ) 

FRp  - o.onol 
f K1  = .001 
tRD  i 0.0 1 
SRS=PV»FYV/FPC 

COMPUTE  LOADINCi  COEFFICIENTS 
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V ( 1 ) = 1 .n/Ar. 

AV ( 1 I = V ( 1 ) - Al  I I ) /D 
Uu  10  I -"2  »\ 

V 1 1 1 =V(  I -1  ) -Al'(  I -1  ) / Ar. 

in  A‘-'(  I ) = Ai-'.  ( I - 1 ) + y ( I 1 - AL  ( i ) / IJ 

iNiTiALi/f  npi  Lt:>^r  luNs 

no  1 1 I = 1 .N 
A A 1 ( I)  - n . n 

A A 1 ( I ) 0.0 

AK2''AX  1 I ) =n.o 
tMAX ( : ) =0*  0 
As  l \',Ax  ( I ) -0  • 0 
AK10X( 1 )=0.0 

DhLUK  ; ) - o.n 
11  oAl  TA ( 1 ) =n.o 
PA]  = o.n 
‘■COAX  = o»n 

xDfLlA  = 0 
I.'LLT  ‘O- -0.0 
FC(*1)  - n.np 
L EC -0 . ^00  1 

cO  : I ( EC  ( \ ) . 'I . E'j*- . 1 ) : .001 

EC  ( iM  =LC  ( : Mj-s 
?3  UEB  = .001 

EB(r.)  = OEB-EsC.) 

CALC  -’LAIL  AXI’.l  IOiXCl  A.'.n  ■IJi'mEAT  AI  beam  Cl  .1E^LI.'.L 
21  call  FOKCEC..!  s!M  O.BC.;  .PA) 

AEOr  15,-.  S { 1V(  ‘.  ) +")EL  T A I \ ) 1 

APPLUC  ■■COM  !.■  KLS.oll^.C)  ■Aor^t  ■,  1 

I f ( .iflO  ((  R''H  0 ) -Al  ■•)/.<•■(•.))  .L  I .1  iXO  ) c,o  Tl'  30 
I F ( ‘2  o ( \ ) . ’ • A E '•! ) o 0 Tv’  22 
EBi  M ^ EB  ( ’J ) 0>Lu 

00  TO  21 

2 2 EB(  0 ) - EE  r;  ) - . ’■J-  Oi  B 

0EH  = r)f  !•/<.  . 
f-2  TO  21 
X''  CO'UIN  'f 

•.O^A  HAVE  PA  AT  CLMt^E:’)!  --  f I *U)  T’A<AME1LRS  AT  uTmLR  NuULS 

'.)  = %_ : 

1 O ‘.00  i-l.v 
J = N-  I 

AE'3=PA*(AM(J)+Df  L TA  IJ  ) ) 

EC2^FC ( Jtl  ) 

IE  ( A<  ,GT  . 1 .0  )f.;)  TO  •.‘‘0 
I E ( I . t (J.  1 . A'.IIj.PA1  ,(,T  . PA  ) LC2  - LCMAX 
APo  call  e c a,'t!jx  ( j ,t  l2  .pa  . AL  M ,F  C ( J ) ,‘  r ( i\  ) ) 
bnn  co.\tinue 

COf-’PUTF  DcFLFC  T I OHS  FkuM  Ntw  PHIS 


I 
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KnCLiA  = \DtLTA-fi 

[J  tnf  t A U 1 = AH  1 I “ l',i  1 i 1 ) / I 3 . »i) ) ^ AL  ( 2 ) » I . Pii  n i ) < I It  I ( ^ ) / / 1 u . »Li  I 
THl  T Al  1 ) I A ( 1 ) 

i)-'  AO  I-2»N 

0 IiU,  I A I 1 > - ' AH  1)  •>  i ' Pm  I t ; I 1 1 ,1 1 1 1 - 1 ) ) + Al  I 1 + 1 ) « H . ‘ t M i l 1 ) + P|1 1 > 1 + i ) ) 

1 ) / I s . »D ) 

*-r  THlT  A ( 1 ) ^ IrtL  T A ( 1 ) *-0THM  A ( i 1 

; f ( rc  ( \ ) .G  r . "u  1 1 Hf  I A I ; ' = I ..t  r a i 1 1 + 1 Ptt  1 1 h -^'h  i ( \-i ) ) / a .o 
Of  L T A ( 1 H TftL  1 M 1 : * AL  ( 1 ) /n 
DO  Al  I =2.!. 

THi.  T A ( ! ) = rm  T A ( I -1  ) -r  IMF  I A t [ - 1 ) 

A 1 OLl  ^ A ( I ) = :■(  L A ( I - 1 ) + r M;  T A ( 1)  *AL  ( N / O 

COMPARt  Mw  Dlt-TA  ro  PRt  ; HGG-  01  L T A 

1 r ( A H > 1 1 D L L r A ( N ) - l) L L I A .‘-J ) / 0 L L T A ( N ) ) , L c.  • L l\  G 1 Go  T 0 u 0 


I 1 ( kllPL  TA 
D^  -01  I 

.GF.6)  GO 

= 1 .i\ 

TO  b02 

60  1 

DFLTP ( I ) 
DEL  tan  = 
(.0  TO  ?6 

= DEL  r A ( I 1 
DELTA ( N ) 

0? 

OO  603  I 

= 1 .N 

T-'ELT  A ( I ) 

= (DELTAll 

l+DELTBI 1 ) ) /2 

60  3 

DFltB 1 I) 

= DEl  TA ( I 1 

1 

,-)Fl  TA\  = '^FLTA  ( N) 
k.DlLTA  = 0 
GO  TO  23 

bOO  01  LT  AN-OtL  I A ( "i  1 

D'.  61  1 1 = 1 .% 

I M 1 C ( I ) . I L . L A X ( I ) ) bU  TO  611 
“^MAX  ( n = LC  ( I ) 

A<1 vA  X ( I 1 = A<  1 ( I) 

AL?va,x(  11  = AK,?  ( 1) 

611  COGTIMLF 
PA=PA*CP 
DO  601  I = 1 .N 

FSl ( I ) = AMAXl IfGl ( I ) ,fS(  I ) I 

ESI  ( I ) : AVAX 1 ( LSI  ( I ) . t S ( i ) 1 
AiLlr^X  ( : ) = AMAX  1 ( Ak  1 0 ( 1 1 ,Ak  1 OX  I i I ) 

R ( I ) = R M { 1 I * C M 

fCTl  ( ! ) = (tCI  1 )»(Akt*(  1 l-O.OH)  )/AKA(  I I 
f CT2 ( I)  = ( lC ( I H ( A<a ( I ) -0. 20  ) ) /ALA { I ) 

OELTDl  I l=rfLTA(  I )*|J 
A ( j ) =v  ( 1 I » PA 

SRC ( I 1 =F5C0M  CrON(  I ) tAk A ( I ) I »CP 

601  SRT  ( I 1 =SR'C  ( I I 6KS«  rP 

DElIIED  OUT  PRINT  STATEMENT  DEC  13.1V73 

PRINT  calculated  BtAv  RESPONSL 

«RIt  E I 0 » 3 M 1 .KVA  1 ) .PA  .LEL  I o I 1 ) . f ME  1 A I I ) . E C ( I ) . LH  1)  • oA  l I I . oK  I I I I . 
1 Sf’C  ( I I . Ak  A I i ) , Ak  10  I 1 ! . 1 = 1 .N  ) 

PRINT  6C2 

602  f 0RM\T  ( ijt  ,/  ) 
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n n 


PA=PA/CP 

!i  (PAl-PAl  6 TA  . 6rw;  , ftp  5 

t^^C'-'AX  = LC(N) 

PAl  = PA 

1 F ( rc  ( ,\l  ) . '-t.  .00^^  ) GO  To  6Pj 
GO  TO  20 

CALL  PL  T5-  ( OK  , I AP  .An  .FPL  ) 

CALL  PLOT  10.  . . »-3 ] 

It  (XXX)  ttOO . do  1 . »02 
CALL  PLOT  ( 0 . , 0 . .-.0  ] 

■FTOP  ICl 
F \0 

FUBRO'JT.NE  L C ANPs  ( I . f C 1 . P A • AL'  .LC.XXX) 

FUBKCOT  1 \L  tCA,\OK  CAlCOLATLS  A COTJCKlTL  STPA1.\  Af.O  l'IOTA.'.ol  TO 
T'TL  MLoTPAL  AkTS  TiiAT  oATloT  It-O  T HL  APPLlo^;  Lu'.i^o 

co”vo\/B/TPi  ,lr:x 
DlC-AoiajI  (LCi/'2.  »*00l  1 
tC=FCi 
■.  = 1 

call  FKAP(L^  .)*A.  1 . ’“"’.XXX) 

< = <■•■  I 

IF  ( ABS(  ( AlX-RV) /Atv)  .Lf  .tLpR)(.L  TO  10 
I r ( <.r,T  . Tip  )GO  ' : 10 
I F ( R-''.LT  . Af-' ) G' ' r ■'  P 
Fl  = F C-PFC 

IFire. GT.n.O)  •_  I-  1 
''FC  ^ "ET/2. 

re  = rrc 

GO  TO  1 
FC=FC+.  ’5»Pt  r 
I F ( r C . oT . X XX ) oO  I i 11 
DFC  = PFC/A. 

GO  TO  1 

ore  = .nopi 

FC  = ECi  • .7-,’DtiC 
PFC  = DEC /A. 

C,0  TO  1 

RFTIJR.X 

Er-ID 

GURROurr.t  lOkKiii.*  .tp.  .pa) 

GI/BROU  F I i'(E  FORlL  CAlLjI.AIfLi  LoGCKLTr.  A'.Q  REHAR  Foi\CtS  blALG  A 
.',TRAI\  DI  STR  :loT  1 


G(jV-‘UP/A/AK3.Aj.’_u»^CK.rPC>t  Vj»t  » .Ftp.'  fV.c.r>'"'.''oH.UOIi.L_>tiR»Aor(H. 
IBSHO  ,P  ,0,  PP  , fXC.A  . Ak  (S  ,A<6  . XPC 

COM'AON/C/RM  ( 20  1 G i 20  1 .AKA  I 20  * . AK  10  ( 20  > • Pfl  120'  .CCoiN  ( ..' J ' . f i ( Avj  ) . 

(Fsi(2o)«r-:(20i  •blagx.  t 20  > . ak  icx  1 20 ) . egpp  ( 20 ) 

enw-xon/F  /AKK  1 ( 20  ) •Ai->  2 ( 20  1 
IF  ( PC  +Fl  . • 0.0.0  ) iGO  TO  30 
AK  A ( I ) = AK()»LC/  ( EC  *EB  ) 

ECRR-ECk 

PHI ( I ) =LC/ Aka ( 1 ) 

ESP  = EC*(Aka(I)-Ak.8)/Aka(I) 
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M n o 


T 


t S(  I ) = EC*  ( 1.-AK,h(  I ) 1 / AH...  ( n 
(■■.O  TO  36 
30  CSP=EC 

r S ( n = - f (' 

AKA ( I ) = 1 on. 

36  C ALL  COreKiT  I tC  . A.,  4 ( I ) ,A6?  , : ! 

S c I ) = S r E E L ( I j I I 1 . K T , i S M » A ‘ . f ' . 0 11 » L ‘■1  ,‘i . E j I i I ) . F S 1 ( i ) ) 

Ak^  = '■r;(  J ) »p/frpc 

AK^o  = STEt_L(L6^''.F'T'o,L,,MP*-^}H^  .c^  tip.;_k*  .0.0* 

IF  ( f AO.  ro  . I .,  ) A'*  7P  = Af  70  XP  ( ( Eo-  1 ’ o : ; . 001  1 

IF  I F B . - L . O . 0 ) &>^  T ! 0 

I F’  ( OH.:,  T . : CR  ) r,*,  l i 

AkR=ECR“A-.A  . I 1 /EC 

AK  1 ■'  ( ! ) - ( Af.o-Al,*  I I ) _Ao  » ! s[ 

I f ( AK  1 0 i I ) . 'oL  . A K i 0 X I I ) ) Go  r u L 0 

\ X JO  - •,<  H - ( A k 10X1  I ) - A L 1 0 ( I ) i / D 

F CRR  = F CK»Ax'^<)/A6  A 
A k,‘)  = X/, ‘V 

I F ( A«  O.  ' : . 0. 0 ) AKR=  0.0 

GO  TO  EO 

10  A*kA  = o*n 

A K I o ( I)  = 0 . C 
GO  TO  20 

11  AK9= Ax6-Ak A ( n 
AKlO  ( I ) •-  0. 0 

70  E=E0RR 

I f 1 FB.LL.LCR  I L = tB 
C ON ( 1 I = A< 1 « A< 3 »Aka ( I ) 

TC0F)  = E»E*'C0N'^A>-9/2  . 

PA  = rC0Ml  I ! +Ak  n>-AK  7-TC. jN 

R'k'i  : l=CC0N(  I )»(AK6-XPC-Ak2*A.s-I  I ) )+Ak/F-'«  IAK6-XRl-AkU)+Ak/«  IaHL-AKO 
1 ♦ 1 . ) F-  ICON*  ( XPC- AK6  + AK'.  ( I)  +2  . « Ak'y/  f . ) 

Ak,k  1 ( I ) = A^  1 
Ak<2 ( I ) =AK2 
ESPP ( I ) =EaP 
RETURN 
END 

S'JBROUT  INE  FKAPIEC.PA.I  . K.V  1 . XXX  ) 

SboRO'jTINE  Fk'.P  CALCULATES  Tul  DISTANCE  TO  THE  NLLTKAL  AXIS  Fo.x 
A GIVEN  CONCRETE  STRAIN  TO  sAIISFv  APPLIED  loAdS 

CO^^'XQN/B/FRR  . I CL 

C0M‘xCN/C/R'-'  (20)  »ES(20)  .AM.(20l  .AKl0(20)  .Rh!  (20'  .CC0N(2u)  .FS(<-D)  * 

IF  SI  ( 20  ) .LS  1 ( 20  ) .i  BFi  ( 20  ) 

CALL  FOROPf  I .EC  .-EC  .DF'AX  ) 

I F ( PA .GE .PMAX ) GC  TO  2 
I COUNT  = -1 
DEB  = .001 
EB  = EB3(  I + n 
1 CAI  L FORCE( I .EC .EB.PAl I 

I F ( AhS ( ( PA-PAl  I /PA ) .LE .ERR  ) GO  TO  10 
:r(7A.LT.PAl)  GO  TO  ‘y 
I COUNT  = 1 

EB  = rB-0i:Fi 

IF  (FC-t-EB.uT.0.0)  GO  TO  1 
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n r»  o ^ 


EB  =PEB-EC 
U'J  TO  1 

9  1 M I COUNT  ) 1 1 . 1 1 . 12 

11  EB  = Lb  [Tl B 
GO  TO  1 

12  EP  = ( B + 0.  »HEB 
,TEP  = nrP/9. 

GO  T&  1 

2 NMI  - n.n 
hiEIUT^N 

10  Kf-il  = KI-111) 

EBP  ( I ) --EB 
RETURN 
ENO 

EuNcTION  blLLL(t_..>«f  i»L.jrl»AoU*ujil*‘Mi»Eb'J«fO*jl 

fUNTTION  STr_t  L i-ALCUEAlL;-  A .STEEL  j^RLS.s  GIvtN  A lIlLl  LTkAI.'* 

F Y = FY/EMS 
F SA  r-Af<S  ( rS  ) 

I  t ( E50.  [ Q.  0.0  ) ^>0  TO  2 
U ( : S.GE.f  SU ) GO  TO  2 
STEEL  = ESO- ( LSO-ES  I ■<  ERiS 
RE  TURN 

2  I f ( ESA .LE.EY  ) GO  TO  10 
lElfiSA.LE.ESMlLi'J  TO  11 

S TLl  A = E Y ♦ A6tu>  ( t Sh-LSU  ) -BSH  » 1 lS  A~lSM  1 * « 2 
GT)  TO  12 

10  .STEE  L = ES»E 
RETURN 

11  RTELArEY 

12  I L ( t S . L T . 0 • 0 ) GO  TO  n 
1 steel  ^ E T E ^ A 

RETURN 

13  STEf  L = -‘^  T EL  A 
RETURN 
END 

EUMCT ION  1 VALUE  I E .LO.Lu,  1 ) 

FUNCTION  value  CALCULATES  A vAlUE  OE  LONCKETL  SIKLSS  u 1 V E N A 
r CONCRF TE  STRAIN 

C 

E'A.  ( rO.EU  ) = 2.  »E0  / ( 2 . • (1  U«  »2-E0«  *2  ) ) 

F2'' - ( 0. 3 + 1 0»FM  ( ( O.E  J ) ) / f K I EO.LU  ) 
iru  .GI.E20)GO  TO  3 
irc'.GT.roGO  10  2 
F VAlUE  = 2 *E /F  0- ( E /EO) > « 2 
GO  TO  10 

2 F' VALLir  = 1 .0-  ( E-EO  ) «F  Y{  EO.EU  ) 

GO  TO  10 

3 fvali;f=o.2 

10  RETURN 
FNO 

.SUElROUT  1 NE  CONCRT  ( A , AK  A . A»f  1 . AlC  2 , 1 ) 


R(ACT  = :.0 
F = A 

I  f ( ' .r,T  .F.'^^X  ( I ) ) A..  lo  1 

F CVAX  -f  VAL'-.'t  ( ! “A  X ( I ) .KSE  , - I ) 

FC  = I f ( EMaX  ( 1 ) - F ) • F L.’ILOD 

I r ( .1  • .0.0  ) Gl.  TC 

RFArT  ^ fC/EC^’AX 
L ^ E“AX(I) 

I F ( ',f;a  . at  . AK  a ) GO  ’ . 20 

AX.l  - RFAv  T «A«:  r-'.  AX  ( I ) 

AF  2-  At  V ( : ) 
pr  T'  lOv 
2 A X 1 - 0 . A 
A t 2 = 0 . ’0 
RFTUPri 

1 I F ( axa/Ak<.  .1  ’ . 1 .0  IC.O  Tt'  2 0 

IF  ( t.  . G i . L 2 0 ) o o TO 

II  ( t-  .r,  I .r  U 1 GO  TO  2 
A K 1 = F A 1 U ) / E 

A>  2=  1 .0-1  HU  (t.  I /L 
K’ETliP.N 

2 AKI =r ARFA2 ( t ) /[ 

AX?  = l .G-F  FF>,T2(L  ) /I 
Of  Tl  iRF! 

3 A<1  = F ART  Ail t ) /i 
AX?r  1 ,0-f‘HU  il  t ) /f 
RF  Tl  iRN 

n ALOhA= 1 . O-AXft/AXA 
f p = p«AL(  HA 
I F I r .AT . f?0 I GO  TO  5 
ifif.gi.f.o)  go  T^  a 
AM  - ( FAl  ( L ) -(  A1  I LR  ) I / t 

Ak2  = 1 .0  - II  A 1 ( L ) »F  LiJl  ( t ) - F A 1 ( LFU  »l  Lol  I lH  ) ) / I ( F Ai  I t ) - F A1  ( LF'  ) ) *L  ) 
liO  TO  loo 

A I F I LP.GTji!  0 ) GO  TO  6 

AMrIFAOEAPILI-FAllEPI  )/£ 


GUHRCtllINF  CONLKI  CALGOLAILj  OuMi,KtlL  bIKLS.S  t5LO^^  RA«Ai-iLTLkP 
GIVEN  A ^.uM-'KLTt  GTIMIM  ANi;  Ol.STANOt  TO  ThL  NlUTKaL  mXIG 

t llMAlG'l  / A/AF.  i.L,0.L^,  ,tCR»''PC.Li''.,j»FY.r  FF’.I  YV»Lbii»/'oit»UoH.LonF^*«oriH, 
IHSHP  ,P  ,D  ,RF'  ,E  MCON,  Ak.6  , AK8  . XF'C 
COVMON/LVt  RAX  ( 20  ' . A<  H-‘AX  ( 20  I . \it2‘AAX  I 20  ' 

F:1(rO.EU)=2.»lG/  I3.*ll  u*»?-L0»*2)  ) 

F A1  (I  )-|  x.<LO*I  «»2-t  U*»lo*«2I 

FFAMF  )r(B.«!  »CO-l.»f»»2l/ll2.»t0-4.xf) 

F A2  IF'  1 = t -0 . 8 »T'’-‘  ( t 0 .tU  1 » * 2 

r Hi  2 n ) = ( ^ ‘ - 2 . " F I F t . F 0 ) “ E « X 2 ) / I 6 . - 3 . H M ( F.  0 . 1 U ) * L ) 

' A 'U  £ ) = 0 . 2 ♦ ( 

F F P X ( r ) = o . 5 » E. 

F ARf  A2  t E ) A 1 l EO ) tFA? I F -tO) 

FLBT2lEI  = lFAl(LiJ'»rtullLO)»-F  A«  (L-LJ)'(Lj  + F-  |(tlL-u>^H  J/FAktA2'LJ 
F AR  t A 0(  E I --  F A 1 I L 0 ) -t  l A 2 U 2 0 - F G I < I A il  E - .2  0 ! 

F t HT  i ( E H I F Aj  ( Eg  I H t ••  1 ( I G ) + F A^  ( l20-‘-G  ) « ( L G*  I Lb2  I L^O-t  G ) Ht  Ai  Ih-^ 

1 ) M r 20  t t t f.  :•  ( L-I  20  ) ) I .'F  ARE  A J I T ) 

EUNL02.--2. 0/20 

f.  2 n H 0 . 8 * E 0 » F '•!  I L 0 . 1 u ) I / 1 I L G . L U ) 


Ak?  = : .0  - ( =■  ^ r-^;  1 L ) 'r  . t-i  T2  ( L ) - f • i I .J  ,<■.£;<  1 1 lii  ) ) / 1 I t mKLA2  I L ) -(-  Ai  ; n 

1 »r ) 

'•  if  20)  'iO  TO  i 

A),  I - ( F A 'U.A2  1 (*)>))  / l: 

,'OAfTOr  ( i F ; I,f  ',2  ( r 1 - --  A.Rf-  -,2  ( r i'  I ) ' L'  ) 
a:  AA2  = F'AiU.  A2  ( L ) « r EHT2  ( r ) 

AF'AP.2-f  A -:r.  A2  ( EP  ) f F'D^  2 t EP  ) 

AA2-  1 .0-(  APAA2-AFAP2  ) //«P/'!.2 
CIO  TL'  ICO 

A IE  .1  20)CjO  Tt;  7 .' 

1 ( ( F P • w T , EC  ) 2)  j T 'v/  3 

AKl  --  ( F APE  A’  ( L ) -i  . 1 i LP  ) ) / L 

A ,■?  7 r ( ( E A ^ A 2 t L ) - f /■■  1 ( L P I I U.  I 

A !A2  " 1 • 0 - ( f*  Al^*  _.A  ^ ( L ) - f i..  t?  1 7 ( t ^ 1 i r ) ^ L L’  1 i L ►-  ) ) / 'A  t-'/.  . ^ 

CO  ’■0  100 

E A'-  1 = ( F APE  A A ( 1.  ) - f ►'PL  A2  { 1 E ) I / L 
APAp  = ( ( , a.K  /•.?(:.)-►  ^L7  2 i lP  ) I L 1 

ArAAt>2  = ‘ APlA3  ( !-  ) ►)  foil  3 ( L ) 

..f'A'?iS2  = F'^  Pl  A2  I E ) • I ti-  T 2 i Lf  ) 

7 .'<.2  = : . 0 - I APAAB2--'.PAiia2  ) / apaK 
O'  - T j 1 1 0 

7 AF.l-^  ( F APLA3  ( L ) -FAP.LA3  ( tP  ) ) /i 

APAK  1 = ( ( f AfetAi  ( ci  ) -(  A.vL  A M i P ) ) n ; 

APAAC7^FAPLA3  ( t.  ) «FEl-  T •'  ( t ) 

A PA  r.C  2 = F A ,■!  L A A ( E P I » r F P T 3 ( L F ) 

AK.2  = 1 . A - ( APA  AC2-APAEC2  ) / Al'AK  1 
oo  TO  inn 
V PKI'n  25 

F'oPoat  ( ///lnx*^f  P Lxcrto'  tPo*) 

STOP 
n A)  1 ■: 

Pt  TUPN 


AKl «KEAC  T 


ABBREVIATIONS,  ACRONYMS,  AND  SYMBOLS 


A , , B . 
sh  sh 


"unload 


area 

factors  relating  to  lateral  load  distribution  in  beams 

area  of  tensile  reinforcement 

area  of  compression  rei nforcement 

constants  relating  to  steel  stress-strain  curve 

compressive  concrete  force 

compressive  reinforcement  force 

modulus  of  elasticity  of  concrete 

modulus  of  elasticity  of  steel 

unloading  modulus  of  elasticity  of  concrete 

lateral  load 

axial-to-lateral-load  ratio  (P/F) 
length  of  beam 
moment  capacity  of  section 
moment  coefficient  at  node  i 
resisting  moment 
resisting  moment  at  node  i 

factor  used  in  expression  for  moment  capacity  of  beams; 
number  of  nodes 

axial  load 

tensile  force  in  concrete 

tensile  force  in  reinforcement 

shear  coefficient  in  segment  i 

depth  of  concrete  tension  stress-block 

distance  from  bottom  of  beam  to  plastic  centroid 

depth  of  rectangular  concrete  stress-block 

shear  span 

shear-span- to-beam-depth  ratio 
width  of  cross  section 

distance  from  compressive  face  of  member  to  neutral  axis 

distance  from  compressive  face  of  member  to  centroid  of 

reinforcing  steel 

compressive  stress  in  concrete 
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ABBREVIATIONS,  ACRONYMS,  AND  SYMBOLS  (Cont'd) 


compressive  strength  of  concrete  test  cylinder 

flexural  strength  of  concrete  in  compression 

f^^  cube  strength 

f . forces  at  node  i 

f stress  in  tensile  reinforcement 

s 

f^  stress  in  compression  reinforcement 

f ultimate  steel  stress 

su 

f^  yield  strength  of  reinforcement 

k ratio  of  distance  from  compression  face  to  neutral  axis 

to  effective  depth  (d) 

ki,kj  stress-block  constants 

ks  factor  relating  strength  of  concrete  in  beam  to  concrete 

cyl i nder  strength 

m slope  of  descending  portion  of  concrete  stress-strain 

curve  used  for  analytical  model 

bending  moment  at  node  i 

n modular  ratio 

p rei nforcement  ratio 

p critical  reinforcement 

^cr 

t beam  depth 

v^  shear  force  in  segment  i 

6^  deflection  of  node  i 

c strain 

c.  strain  at  bottom  of  beam 

b 

c concrete  strain 

c 

L . concrete  strain  at  node  i 

Cl 

concrete  strain  at  cracking;  limiting  concrete  tensile  strain 

£ concrete  strain  at  maximum  concrete  stress 

0 

strain  in  tension  reinforcement 

. ^ strain  in  compression  reinforcement 

■ strain  at  commencement  of  strain  hardening  in  steel 

t steel  strain  at  maximum  stress 

su 

ultimate  concrete  strain;  strain  at  maximum  ki 
£?o  concrete  strain  at  stress  of  20  percent  of  maximum 


ABBREVIATIONS,  ACRONYMS,  AND  SYMBOLS  (Concl'd 


plasticity  ratio 
rotation  of  node  i 
ducti I i ty  ratio 
collapse  ductility  ratio 
curvature  at  node  i 
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